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Abstract 



The production of neutrons carrying at least 20% of the proton beam energy {x L > 0.2) 
in e + p collisions has been studied with the ZEUS detector at HERA for a wide range 
of Q 2 , the photon virtuality, from photoproduction to deep inelastic scattering. The 
neutron-tagged cross section, ep — > e'Xn, is measured relative to the inclusive cross sec- 
tion, ep — > e'X, thereby reducing the systematic uncertainties. For x L > 0.3, the rate of 
neutrons in photoproduction is about half of that measured in hadroproduction, which 
constitutes a clear breaking of factorisation. There is about a 20% rise in the neutron 
rate between photoproduction and deep inelastic scattering, which may be attributed to 
absorptive rescattering in the jp system, or 0.64 < x L < 0.82, the rate of neutrons is 
almost independent of the Bjorken scaling variable x and Q 2 . However, at lower and 
higher x L values, there is a clear but weak dependence on these variables, thus demon- 
strating the breaking of limiting fragmentation. The neutron-tagged structure function, 
F\ n(3) (x, Q 2 ,x L ), rises at low values of a; in a way similar to that of the inclusive F2(x, Q 2 ) 
of the proton. The total 77r cross section and the structure function of the pion, F^ix^, Q 2 ) 
where x n = x/(l — x L ), have been determined using a one-pion-exchange model, up to 
uncertainties in the normalisation due to the poorly understood pion flux. At fixed Q 2 , 
F£ has approximately the same x dependence as F<i of the proton. 
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1 Introduction 



Leading baryon production in a process with a hard scale provides a probe of the relationship 
between the quantum chromodynamics (QCD) of quarks and gluons and the strong interaction 
of hadrons. The leading baryons are produced with small transverse momentum (p T ), guaran- 
teeing the presence of a soft process with its related long-range correlations. The hard scale can 
come from large photon virtuality {Q 2 ) in deep inelastic scattering (DIS), or high transverse 
energy (E T ) in photoproduction. The observation of events with neutrons or protons carrying 
a large fraction (x L ) of the incident proton beam energy in positron-proton (e + p) scattering at 
HERA has led to renewed interest in the QCD evolution and factorisation properties of 
proton fragmentation to leading baryons in deep inelastic scattering [|7|-[rT| . 



A non-perturbative approach to the strong interaction is needed to calculate the production 
rates of leading baryons. Some theoretical models retain the QCD building blocks, quarks 
and gluons, as fundamental entities but add non-perturbative elements, such as soft color 



interactions [|L2 |. Another approach interprets the production of forward baryons in terms of 
the exchange of virtual particles [|T3|-f23|j , which accounts for charge-exchange processes (p — > n) 
in hadronic interactions |24H33|. In leading neutron production, the large diffractive peak in the 
cross section, which is observed in leading proton production j|,|5]], is absent. Although leading 
neutrons can be produced indirectly through the production and decay of heavy baryons or 
through the fragmentation of diffractively dissociating protons, direct production dominates || . 
In the language of particle exchange models, isovector exchange (tt, p and a 2 ) is required 
for direct neutron production, and isoscalar exchange (/ and uj) is absent. The production 
mechanism of leading neutrons is therefore simpler than that of protons and, since the cross 
section for 7r exchange is much larger than that for p or a 2 exchange, it is predominantly the 
structure of the pion that is probed by the hard process 



Comparisons between cross sections for the production of particles in the fragmentation region 
of a target nucleon provide tests of the concepts of vertex factorisation and limiting fragmen- 
tation. Factorisation tests involve comparing semi-inclusive rates, normalised to the respective 
total cross section, to study whether particle production from a given target is independent 
of the nature of the incident projectile. Thus, for example, data from jp interactions may be 
compared with those from pp collisions. The limiting- fragmentation hypothesis ITDfl states that, 



in the high-energy limit, the momentum distribution of fragments becomes independent of the 
incident projectile energy. This has been verified in hadronic interactions |34j| . It has not been 



extensively studied for produced baryons in electroproduction. In this case, the fragmentation 
spectra should be studied as a function of Q 2 and the photon-proton centre-of-mass energy, 



W 11 



In electroproduction, the size of the incident projectile (the virtual photon) can be varied. At 
Q 2 = 0, the photon has a typical hadronic size; as Q 2 increases, the photon size decreases. 
Hence, absorptive rescattering of the produced baryon is expected to decrease and become 
independent of Q 2 at high Q 2 |35|,|36|]. Such effects, if x L or p T dependent, lead to an apparent 
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failure of factorisation and limiting fragmentation for low values of Q 

This paper reports studies of the production of leading neutrons using the ZEUS detector and 
its forward neutron calorimeter (FNC). Leading neutron production can be measured precisely 
by determining the neutron-tagged cross sections relative to the inclusive ep cross section, 
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which has been measured with high accuracy |J7|. The study of ratios reduces the systematic 
uncertainties arising from both theoretical and experimental sources. 

The paper is organized as follows: the kinematics of leading neutron production is defined in 
Section 2; the one-pion-exchange model is reviewed in Section 3; the detector, the reconstruction 
of the kinematic variables and the experimental conditions are described in Section 4; the event 
selection is discussed in Section 5; in Section 6, the advantages of analysing the data using the 
ratio of the neutron-tagged to the inclusive cross section are outlined; systematic uncertainties 
are discussed in Section 7. The general features of the neutron energy spectra and the ratio of 
leading-neutron to total cross sections are discussed in Section 8, where the factorisation and 
limiting-fragmentation properties of the data are investigated and an estimate of absorptive 
effects is made. The neutron-tagged structure function is discussed in Section 9. The data are 
analysed in Section 10 in terms of a one-pion-exchange ansatz and the pion structure function, 
F£, is extracted up to uncertainties in the normalisation due to the lack of knowledge of the 
pion flux; theoretical models are compared to F£. Section 11 summarizes the results. 



2 Kinematics and cross sections 



Figure p] shows a Mandelstam diagram of semi-inclusive leading neutron production in ep col- 
lisions. Four kinematic variables are needed to describe the interaction ep — ► e'Xn. They are 
defined in terms of the four-momenta k, P, k' and N, respectively, of the incoming and outgoing 
particles e,p,e' and n. Two variables are chosen from the Lorentz invariants used in inclusive 
studies: Q 2 = —q 2 = —{k — k') 2 , the virtuality of the exchanged photon; x = Q 2 /(2P ■ q), the 
Bjorken scaling variable; y — q-P/{k-P) c± Q 2 /(sx), the inelasticity; and W 2 = (P + k — k') 2 = 
m 2 + Q 2 (l — x)/x, the squared mass of the produced hadronic system. In these equations, m p is 
the mass of the proton and y/s = 300 GeV is the ep centre-of-momentum-system (cms) energy. 

Two further variables are required to describe the leading neutron. They can be chosen as the 
laboratory production angle of the neutron, 9 n , and the energy fraction carried by the produced 
neutron, 

N ■ k E n 

L p-k~ v 

where E p is the proton beam energy and E n is the neutron energy. The transverse momentum 
of the neutron is given by p T ~ x L E p 6 n . In terms of these variables, the squared four-momentum 
transfer from the target proton is 

t = (P- N ) 2 ^-^-t , 
where t is the minimum kinematically allowed value of \t\ for a given x L 



, / \ (1 x l) ( 2 2\ (1 X l) 2 
t (x L ) = (m n - x L m ~ m , 

x, ' v., 1 



and m„ is the mass of the neutron. 
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The differential cross section for inclusive ep — > e'X scattering is written in terms of the proton 
structure function, F 2 (x, Q 2 ), as 

J2 ep^e'X 

^-^ =m{x , Q ^ 1 + 6), (1) 

where 



47ra 2 / y T 

and 5 is the correction to the Born cross section for photon radiation, Z° exchange and the 
longitudinal structure function, F^. In analogy to this, the differential cross section for leading 
neutron production, ep —>■ e'Xn, is written as 

Ji—ep—^e' Xn 

a ° ^LN(4)/ 



, „ v , , /CF 2 LN(4, (x,Q 2 ,x L ,p T )(l + 5 LN ). (2) 
dx dQ z dx L dp T 

The term 5 LN is analogous to 5, but for the case of leading neutron production. 

In the results presented here, 6 n (and hence p\ and t) is not measured. Integrating Eq. (|2]) 
up to the maximum experimentally accessible angle, #™ ax , corresponding to a p™ ax which varies 
with x L , leads to 

d 3 a ep^e'Xn ^ ^ep^e'Xn 

-dp T 



dx dQ 2 dx L Jo dx dQ 2 dx L dp 7 

= /CF 2 LN(3) (x,g 2 ,x L )(l + 5 LN ) (3) 

where 

max 

F 2 LN(3) (x,Q 2 ,x L ) = / T F^ 4) (x,Q 2 ,x L ,p T )dp T 
Jo 

is the neutron-tagged structure function integrated over the measured range in 8 n . 

It is sometimes more convenient to discuss the results in terms of cross sections rather than 
structure functions. At low x, where W 2 ~ Q 2 /x, the following expressions are valid: 

*fiWQ 2 ) = ^?^(*,Q 2 ) (4) 

da^ p ^ Xn (W, Q 2 ) fPT x d 2 a^ p -* Xn (W, Q 2 ) , 

~~ ' — dp T 



dx L Jo dx L dp T 

= ^F^( Xl Q 2 ,x L ), (5) 

where cr^ ot p is the total virtual photon-proton cross section and d 2 a' 1 * p ~* Xn (W, Q 2 )/dx L dp T is the 
differential cross section for leading neutron production in a virtual photon-proton collision. 



3 One-pion-exchange model 



A successful description of the available data 



on charge-exchange processes (p — > n) 



in hadron-hadron interactions has been obtained using the exchange of virtual particles with 
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the quantum numbers of the tt, p, and a 2 mesons |T3|-|iT| ^Hj, |2"T|] . In such processes, the pion, 
due to its small mass, dominates the p — > n transition amplitude, with its relative contribution 
increasing as \t\ decreases. In contrast to the nucleon, whose partonic structure has been probed 
directly in deep inelastic scattering, the structure of mesons has been studied only indirectly 
using hadron-hadron collisions. For the pion, the data from np interactions yielding prompt 
photons |3l| 55 1, high-mass muon pairs |4"0"|-44"5|] and dijets |f0] in the final state have given 
important information. However, these results are mostly sensitive to the high-x, or valence, 
structure of the pion: the sea-quark and gluon content of the pion at low x have not yet been 
measured. The HERA leading neutron data offer an opportunity to probe this low-x structure 
via the diagram shown in Fig. [j]. In the one-pion-exchange (OPE) approximation, the cross 
section for hadroproduction, hp — > Xn, can be written as 



da 



hp^Xn 

dx L dt 



Ip 



X 



L 3 



(6) 



where the flux factor, f n / p (x L , t), describes the splitting of a proton into a im system, s' = s(l — 
x L ), where s' and s are the squares of the cms energy for the hn and hp systems, respectively. 
Charge-exchange processes at the ISR [27, 28] and Fermilab [|29].[I(]] are well described by such 



an expression with Hl5|, |iq1 



1 2g 



2 

■Kpp 



An An (t - ml) 



1 - x, 



vl-2a! T (t) 



(7) 



where gl pp /(An) = 14.5 is the n°pp coupling constant and a w (t) = a't with a' = 1 GeV -2 . The 
form-factor F(x L ,t) parametrises the distribution of the pion cloud in the proton and accounts 
for final-state rescattering of the neutron. 

For the hadronic charge-exchange experiments, a good description of most of the pn — > Xp 
data 0] is obtained using the Bishari |15|] flux, which corresponds to Eq. ([F[) with F(x L , t) — 1. 
In this case, 

1 2g 2 —t 

f e s(x Lj t) = —- — fl - X, 



An An (t - ml) 2 

The flux-factor / e g can therefore be interpreted as an effective pion flux that takes into account 



processes that occur in hadronic charge-exchange processes, such as absorption (Section |10.1| ) 
non-pion exchange and off-mass-shell effects. 



4 HERA and the ZEUS detector 

During 1995-97, the HERA collider operated at ^Js ~ 300 GeV with positrons of E e = 27.5 GeV 
colliding with protons of E p = 820 GeV. In addition to the 175 ep colliding bunches, 21 unpaired 
bunches of protons or positrons allowed the determination of beam-related backgrounds. 



4.1 The ZEUS detector 

The ZEUS central detector is described in detail elsewhere [|45[| . Here the main components 
used in the present analysis are described briefly. The central tracking detector (CTD) |f46| , 
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positioned in a 1.43 T solenoidal magnetic field, was used to establish an interaction vertex 
with a typical resolution along (transverse to) the beam direction^] of 0.4 (0.1) cm. 

The high-resolution uranium-scintillator calorimeter (CAL) |47]] consists of three parts: the 
forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part is subdivided 
transversely into towers and longitudinally into one electromagnetic section (EMC) and either 
one (in RCAL) or two (in FCAL and BCAL) hadronic sections (HAC). The smallest subdivision 
of the calorimeter is called a cell. The CAL relative energy resolutions, as measured under test- 
beam conditions, are a(E)/E = 0.18/y/E for electrons and a(E)/E = 0.35/a/E for hadrons 
(E in GeV). 

In addition to the CAL, two small electromagnetic sampling calorimeters, situated beside the 
beam pipe, were used to measure the scattered positron. The beam pipe calorimeter (BPC) |48|j 
was located at Z = —2.9 m from the interaction point, and the LUMI calorimeter ||49|| , which 
forms part of the ZEUS system for monitoring the luminosity, was located at Z = —35 m. 
These additional calorimeters had relative energy resolutions of a(E)/E = 0.18/ '-\[E. 



4.2 The forward neutron calorimeter 



A forward neutron calorimeter (FNC) was installed in 1995 |50[ in the HERA tunnel at Z = 
+106 m and at zero degrees from the incoming proton direction. Figure ^| shows the layout of 
the FNC and that of the six stations (S1-S6) of the leading proton spectrometer (LPS). 

The structure of the calorimeter is shown in Fig. |3|. It is a finely segmented, compensating, 
sampling calorimeter with 134 layers of 1.25 cm-thick lead plates as absorber and 0.26 cm- 
thick scintillator plates as the active material. The scintillator is read out on each side with 
wavelength-shifting light guides coupled to photomultiplier tubes (PMTs). It is segmented 
longitudinally into a front section, seven interaction-lengths deep, and a rear section, three 
interaction-lengths deep. The front section is divided vertically into 14 towers, each 5 cm high. 
The relative energy resolution for hadrons, as measured in a test beam, was a/E = 0.65/v^E'- 
As seen in Fig. [|, the FNC completely surrounds the proton beam, which passes through a 
10 xlO cm 2 hole in towers 11 and 12. Three planes of scintillation counters, each 70 x 50 x 2 
cm 3 , are located 70, 78, and 199 cm in front of the calorimeter. These counters completely 
cover the bottom front face of the calorimeter and are used to identify charged particles and so 
reject neutrons that interact in front of the FNC in inactive material such as magnet support 
structures, the beam-pipe wall and the mechanics and supports of the LPS. 

The HERA magnet apertures limit the FNC acceptance to neutrons with production angles 
less than 6 l ,™ ax = 0.8 mrad, that is to transverse momenta p™ ax = E n 6™ x = 0.656 x L GeV. Only 
about one quarter of the azimuth extends to #™ ax , as can be seen from the outline of the aperture 
shown in Fig. The overall acceptance of the FNC, taking account of beam-line geometry, 
inactive material, beam tilt and angular spread, as well as the angular distribution of the 
neutrons, is ~20% at low x L , where the p T range covered is small, but increases monotonically, 



1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the proton 
beam direction, and the X axis pointing left towards the centre of HERA. The coordinate origin is at the nominal 
interaction point. The polar angle 8 is defined with respect to the positive Z-direction. The pseudorapidity, r/, 
is given by — ln(tan |). 
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exceeding 30% at high x L , as seen in Table [I|. The dominant influence on these values is the 
lack of complete azimuthal coverage illustrated in Fig. |j. The kinematic region in t and p\ 
covered by the FNC as a function of x L is shown in Fig. |5](a) and (b), respectively. Although 
the acceptance extends to p 2 T ~ 0.4 GeV 2 , the mean value of p 2 for the accepted data is less 
than 0.05 GeV 2 |]. 

Charged particles produced at small angles and energies ^ 450 GeV are swept away vertically 
by the HERA bending magnets. Higher-energy positively charged particles in the energy range 
from about 450 GeV to about 600 GeV are bent into the top towers above the beam pipe. 
Neutrons are easily distinguished from protons, since the latter deposit most of their energy in 
the top four towers of the FNC. Protons with energies ^ 600 GeV do not exit the beam pipe 
in the vicinity of the FNC. 

Photons, which populate only the lower x L range of the data, are separated from neutrons using 
the energy-weighted root-mean-square vertical width of the showers. After this separation, the 
correction for photon contamination in the final neutron sample is negligible. 

Timing information from the FNC was also used to eliminate backgrounds. One source of back- 
ground is due to the random overlap of energetic neutrons from proton beam-gas interactions 
with genuine DIS events. The rate of FNC energy deposits above 400 GeV was monitored and, 
from comparisons with the bunch-crossing rate, the overlap was estimated to be less than a 
few percent and was ignored. This background was also checked using randomly triggered 
events and events with an energy deposit in the FNC greater than 1000 GeV. 



The calibration and monitoring |5T| of the FNC followed the methods developed for the test 
calorimeters fl|,|5"2|. The gain of each PMT was obtained by scanning the FNC with a 60 Co 



radioactive source. Changes in gain during data taking were monitored using energy deposits 
from interactions of the HERA proton beam with residual gas in the beam pipe. The overall 
energy scale was set from the kinematic end-point of 820 GeV by fitting data from proton beam- 
gas interactions with energy greater than 600 GeV to the shape of the spectrum expected from 
one-pion exchange [15|,|53|]. The accuracy (± 2%) of the energy-scale determination was limited 



by the uncertainty of the shape of the beam-gas spectrum at the kinematic endpoint. 



4.3 Experimental conditions 

The data were taken in three ranges of Q 2 . The first, consisting of photoproduced events 
where the scattered positron was detected in the LUMI calorimeter, is characterised by pho- 
ton virtuality Q 2 < 0.02 GeV 2 ; the second, consisting of intermediate-Q 2 events where the 
scattered positron was detected in the BPC detector, is characterised by photon virtuality 
0.1 < Q 2 < 0.74 GeV 2 ; and the third, originating from deep inelastic scattering, where the 
scattered positron was detected in the CAL, is characterised by photon virtuality Q 2 > 4 
GeV 2 . The three samples are denoted, respectively, as PHP, intermediate-Q 2 and DIS. 

The data for the PHP sample come from four short dedicated runs in 1996, taken with a 
minimum-bias trigger, and correspond to an integrated luminosity of 49 nb" 1 covering the 
full x L range. For x L > 0.5, 2 pb _1 of data were also taken in 1995 with the FNC included 
in the trigger. The other data samples were taken in 1995-97, and correspond to integrated 
luminosities of approximately 14 pb _1 for the intermediate-Q 2 and 41 pb -1 for the DIS samples. 
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4.4 Trigger conditions 



A three-level trigger was used to select the events. The trigger decision relied primarily on 
the energies deposited in the CAL, the BPC or the LUMI and was based on electromagnetic 
energy, total transverse energy and missing transverse momentum. CAL timing cuts were used 
to reject beam-gas interactions and cosmic rays. No FNC requirements were placed on the 
trigger, except for part of the PHP sample in which events were triggered by energy deposits 
in the FNC, LUMI and RCAL. This data sample is useful only for x L > 0.52, where the FNC 
triggering efficiency is better than 50%. The efficiency is 100% for x L > 0.7. 



4.5 Reconstruction of kinematic variables 

The reconstruction of the lepton variables depends on the device in which the scattered positron 
was detected. The LUMI measures only the energy and the horizontal deflection of the scattered 
positron: such particles have scattering angles 9 e < 6 mrad, so that the virtuality of the 
exchanged photon is limited to Q 2 < 0.02 GeV 2 , with a mean value ~ 2 • 10~ 3 GeV 2 . In 
the BPC, both the energy and the position of the scattered positron were measured, and the 



kinematic variables estimated from these measurements |48fl . When the positron was detected 
in the CAL, the double angle (DA) method j54j was used to reconstruct the DIS variables from 
the scattering angles of the positron and the hadronic system. The latter was determined from 
the hadronic energy flow measured in the CAL. 

The impact position of the neutron on the FNC was not measured with sufficient accuracy to 
be used event-by-event. All distributions were obtained by integrating over the angle of the 
produced neutrons up to the maximum observable angle, #™ ax = 0.8 mrad. 



5 Data selection 

The PHP, intermediate-Q 2 and DIS samples were first selected without any FNC requirement. 
The PHP events had to satisfy criteria on the CAL energy deposits or CTD tracks f|9| and 



required a scattered positron in the LUMI detector with 10 < E' e < 18 GeV, where E' e is 
the energy of the scattered positron. The intermediate-Q 2 and DIS events had to satisfy 
35 < 5 < 65 GeV, where 5 = J2(E — Pz) = J2E(1 — cos 9), with the sum running over all 
calorimeter cells; E and 9 are the energy and polar angle, respectively, for a CAL cell. For the 
intermediate-Q 2 data, 5 also includes 2E' e . The events also had to have a reconstructed vertex 
with —40 < Z < 100 cm for the intermediate-Q 2 sample and \Z\ < 40 cm for the DIS events. 
The intermediate-Q 2 events had to have a scattered positron of E' e > 7 GeV within the fiducial 
volume of the BPC Q, while the DIS events had to have a scattered positron of E' e > 10 GeV 
within the fiducial volume of the CAL [E7, 55|. 



The final event selection required, in addition to a good LUMI, BPC or CAL positron as defined 
above, a good FNC neutron candidate satisfying the following conditions: 



• between 164 GeV and 820 GeV (0.2 < x L < 1.0) of deposited energy in the FNC; 
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• energy less than that corresponding to one minimum-ionising particle deposited in the 
counter farthest^ from the front of the FNC, to reject neutrons that shower in front of 
the calorimeter; 

• no signal from the veto counter consistent with a shower from a previous bunch to reject 
pile-up energy deposits; 

• the maximum energy deposit should be in towers 6-9 in the bottom front section of the 
FNC, to reject protons that are bent into the top towers by the vertical bending-magnets 
in the proton beam line; 

• a shower with both energy and vertical width consistent with that of a hadron, to reject 
electromagnetic energy deposits^]. 

The final neutron-tagged sample consisted of 912 PHP, 8168 intermediate-Q 2 and 60102 DIS 
events selected from 56960, 100599 and 642015 inclusive events, respectively. The photoproduc- 
tion background was 1.7% for the DIS sample and 3% for the intermediate-Q 2 sample; since, as 
discussed below, ratios of neutron-tagged and inclusive cross sections are used in this analysis, 
this contamination has a negligible effect. 

The data were divided into x L bins in the range 0.2 < x L < 1.0, as shown in Table [l]. In each 
bin, the acceptance extends from to 0.8 mrad in neutron production angle, representing a p T 
range from to 0.656 x L GeV. 

Different binnings were used for the three data sets. For the PHP data, a single bin was chosen 
with Q 2 < 0.02 GeV 2 and 0.345 < y < 0.636, corresponding to a W range of 176 < W < 240 
GeV with a mean < W > = 209 GeV. For the intermediate-Q 2 data (0.1 < Q 2 < 0.74 GeV 2 ), 



the same 34 bins were used as for the ZEUS inclusive 7*p cross-section measurements ||48|| . For 
the DIS data (Q 2 > 4 GeV 2 ), 25 bins were used, as shown in Appendix I. 



6 The ratio method 

To minimise systematic uncertainties, it is advantageous to use the ratio of cross sections rather 
than the absolute cross sections themselves. The primary measurement, the fraction of events 
with a leading neutron, is experimentally robust: a luminosity measurement is not required; 
uncertainties due both to the finite acceptance of the central calorimeter and to triggering 
inefficiencies are minimised; accurate Monte Carlo modeling of the inclusive hadronic final 
state is not needed; and the effects of radiative corrections are reduced. 

The measured events were binned in Q 2 . Depending on the Q 2 region, the events were further 
binned in x or y. The ratio r(x,Q 2 ,x L ) is defined as 

2 d 3 a ep ^ e ' Xn /dx dQ 2 dx L n ohs (x, Q 2 , x L ) 

r{x : y , x L ) - d2aep ^ e , x/dx dQ2 - A(x L )N ohs (x, Q 2 ) ' 1 j 



2 The farthest counter, 199 cm from the FNC, was chosen to minimise backsplash (albedo) corrections. 

3 This cut is possible since the two-dimensional distribution of energy and shower width shows a clear separa- 
tion between the wide, high-energy hadronic showers and narrower lower-energy deposits arising from photon- 
induced showers. 



s 



where iV bs(a;, Q 2 ) is the number of events observed in a bin of (x, Q 2 ) and n Q b s (a;, Q 2 , x L ) is 
that subset of events with a neutron in an x L bin of width Ax L ; A(x L ) is the acceptance and 
reconstruction efficiency for leading neutrons and is discussed in the next section. 

Inspection of Eqs. ([!]), @ and @ shows that, after appropriate bin centering, 

Fr 3) (x,Q 2 ,x L ) = r( ^ XL) F 2 (x,g 2 )(l + A), (10) 

where 

A ~ 5 - 5 LN - 5 det . 

The additional correction term, 5det, is the difference between the acceptance of the main 
detector for all events and for the subset of events with a leading neutron. The overall correction, 
A, is small and was neglected (see Section |7.3|). 



Equation ( PH| ) can be expressed in terms of cross sections by using Eqs. (Q) and 

dafP^ Xn (W,Q 2 ) r{W,Q 2 ,x L 



dx L Ax, 



-°lt P (W,Q 2 ). (11) 



This equation is also valid for photoproduction (Q 2 =0). It is evident from Eq. (^j) that, within 
an (x, Q 2 ) bin, r is the ratio of the cross section for producing a leading neutron with 9 n < #™ ax 
in a given x L bin to the total cross section. 



7 Corrections and systematic uncertainties 
7.1 Acceptance and corrections 

For the most part, trigger inefficiencies and acceptance uncertainties in the main ZEUS de- 
tector are unimportant for the present analysis since ratios of FNC-tagged events to inclusive 
events are measured. The observed distributions integrated over the (x, Q 2 ) plane are not a 
faithful reproduction the true distributions because of varying trigger configurations and finite 
acceptances. Therefore, when integrated distributions were studied, the events were weighted, 
bin-by-bin, by the inclusive cross section, to produce the correct Q 2 and x (or y) distributions. 
Since the ratio of the neutron-tagged sample to the inclusive sample is not strongly dependent 
on x or Q 2 (see below), this weighting procedure does not significantly change the shape of any 
of the distributions of the ratio. However, the statistical uncertainty is increased since low-Q 2 
events are under-represented in comparison to high-Q 2 events. 

The intermediate-Q 2 and DIS data, triggered by the scattered positron, have a trigger efficiency 
that is independent of the details of the final state. However, the photoproduction data must 
also be corrected for the CAL trigger efficiency, which is dependent on the final state. For 
example, exclusive vector meson production, ep — > eVp, has a low trigger efficiency and no 
neutrons in the final state. As a result, the measured fraction of events tagged with a leading 
neutron for photoproduction must be corrected downwards. A study of each exclusive final 
state and its corresponding trigger efficiency showed that the measured fraction of tagged 
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events in photoproduction should be reduced by (5 ± 2.5)%. This reduction was applied to the 
photoproduction data presented in this paper. 

The acceptance and reconstruction efficiency, A(x L ), of the FNC was obtained using the full 
simulation of the ZEUS detector based on GEANT 3.13 [EfJ. The simulation includes the effect 
of inactive material along the neutron path. The tilt (< 70 /xrad in X and Y, but varying from 
year to year) and divergence (« 70 (100) /xrad in the X (Y) direction) of the proton beam were 
also simulated. 



The relative uncertainty on the absolute energy scale of the FNC is ±2% [^TJ. This introduces 
the normalisation uncertainties shown in Table [I]. 

Because of the poor position resolution of the FNC, the angular distribution of the neutrons 
could not be measured directly. Instead, it was constrained by the distribution of measured 
X positions of the neutrons. The X position of energy deposits in the FNC was determined 
from the ratio of pulse heights in the two readout channels. For the 1995 data taking, the 
corresponding rms resolution was 3.2 cm. The resolution degraded in subsequent years because 
of radiation damage to the scintillator. For fixed x L , the p\ distribution was expressed as 



dN 



dp 



oc e - b{xL)p T. (12) 



The slope b that best represents the data was determined by binning the X-position data in 
x L and comparing it to the Monte Carlo expectation for different b slopes. The results for 
the 1995 DIS data set, given in Table ^| and shown in Fig. |], indicate that the exponential p 2 T 
slope for neutron production rises with increasing x L and may be parameterised by b(x L ) = 
(16.3 x L — 4.25) GeV -2 , as indicated by the solid line (for x L < 0.26, 6 = was used for the 
the acceptance corrections). Also shown in Fig. || as the dashed curve, is the slope parameter 
resulting from the effective OPE flux of Eq. (|8]). The predicted slopes are in reasonable accord 
with the data. 

The polar-angle acceptance of the FNC extends from to 0.8 mrad with an azimuthal accep- 
tance shown in Fig. [|. All of the results presented in this paper have been corrected for this 
azimuthal acceptance. Since p™ ax = 0.656 x L GeV, the sensitivity of the acceptance to the b 
slope decreases rapidly as x L decreases. As a systematic check, the above value of b(x L ) was 
varied by ±2 GeV -2 , which is the approximate uncertainty in b. The uncertainty on the ac- 
ceptance reaches 6% at the highest x L values. The acceptances and their uncertainties, which 
are highly correlated between all of the bins, are summarised in Table [l]. 



7.2 Systematic uncertainties 

The corrections applied to the acceptance lead to normalisation uncertainties. The largest 
corrections and their typical values and uncertainties are: 

• the increased acceptance of the FNC when the LPS is not inserted into the beam, since 
less material is present (—5 ± 1%); 

• the veto cut, which avoids pile- up effects but removes events that should have been 
counted (+7 ±2%); 
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• shower albedo from the face of the FNC, which triggers the veto counters and eliminates 
good events (+1.5 ± 1%); 

• noise in the veto counters associated with the beam (+3 ± 1.5%); 

• veto-counter inefficiency (—1.5 ± 1.5%); 

• the random overlap of beam-gas events with ep interactions (—2.4 ± 0.2%). 

The total normalisation uncertainty for the intermediate-Q 2 and DIS data is ±4%. The cor- 
responding uncertainty for the PHP data is ±5%, which is highly correlated with the DIS 
uncertainty. 

The photoproduction data taken with the FNC trigger were corrected for the trigger inefficiency 
and were used only for x L > 0.52, the region where the trigger efficiency is above 50%. In this 
region, they were normalised to the data collected without the FNC-trigger requirement. This 
leads to an additional normalisation uncertainty for these data of ±4%. 



7.3 Relative corrections 

The relative-correction term, A, of Eq. ( |T0"D accounts for differences in the corrections for the 
neutron-tagged and inclusive samples. Only the difference between the radiative correction 
for the neutron-tagged sample and that for the inclusive sample, (5 — <5 LN , is relevant in this 
analysis^]. The program HECTOR |57], which provides a leading- logarithmic estimate of the 



radiative correction to the Born term in ep scattering, was used to calculate the relative- 
correction term by modeling inclusive neutron production at large x L via one-pion exchange. 
The difference in the neutron-tagged and inclusive correction terms increases with x and x L , 
but remains small throughout most of the kinematic region, reaching (3-4)% at x — 0.06 for 
x L > 0.9 and at x — 0.1 for x L > 0.8. 

The final term in A is due to the correction factor, 8^, which accounts for the differences in 
the central-detector acceptance, migrations and bin-centering for all events and for those events 
with a tagged neutron. 

Since A is smaller than the overall systematic uncertainties of the FNC, it is neglected in the 
present analysis. 



8 General characteristics of the data 



8.1 Neutron energy spectra 

The neutron yields in bins of x L and Q 2 are given in Table |] for the three data samples. 
The data are corrected for the azimuthal acceptance and are integrated over production angles 
< 6 n < 6*™ ax , where #™ ax = 0.8 mrad. As discussed earlier, this fixed angular range corresponds 

4 Possiblc effects due to the longitudinal structure function and Z° exchange are also assumed to be negligible. 
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to a p T range which is x L dependent. The spectra are normalised to the measured number of 
events in the same (x, Q 2 ) bin but without the leading-neutron requirement. Modifying Eq. 
(H) to take corrections into account yields 

nx, v ,x L) = AM ExEq2Ncott ^q2 ) > 

where the subscript indicates that the observed number of events has been corrected bin by bin 
for the effect of the trigger configurations discussed in Section ^. The sum is over the chosen 
range of Q 2 and all measured bins of x in that range. 

Figure |7|(a) shows the x L spectra, as dr/dx L , for the PHP, intermediate-Q 2 and DIS regions. 
Between 5% and 10% of the inclusive events have a leading neutron in the measured kinematic 
range. The ratios are similar in shape for these three kinematic regions, although the fractional 
yield of neutrons slowly increases with Q 2 . The neutron spectra increase slowly with x L (mostly 
due to the increased p T acceptance), reach a broad peak near x L = 0.75, and then rapidly 
decrease to zero at x h — 1. The neutron rate increases by about 20% between the PHP and 
DIS data. 

Figure |7|(b) shows the DIS neutron spectra in three bins of Q 2 . Although, for x L < 0.8, the 
high-Q 2 data tend to lie above the low-Q 2 data, the increase in neutron yield with Q 2 is much 
less pronounced than that seen in Fig. 0(a). In the region 0.64 < x L < 0.82, the neutron yield 
in DIS is approximately independent of Q 2 . For the highest-x L points, the high-Q 2 data lie 
below the low-Q 2 data. 



8.2 Determination of a 1P ^ Xn 

Figure §(a) shows the differential cross section, da hN /dx L , for the photoproduction reaction 
7P — > Xn at a mean jp cms energy (W) = 207 GeV. This result is obtained using Eq. (|TTD 
with a t 7 p t = 174 ± 1 ± 13 jib at (W) = 209 GeV f>8|. Integrating r over the range 0.2 < x L < 1.0 



gives 

r = 5.72 ± 0.12(stai.) ± 0Al{syst.) %, 
corresponding to an observed cross section in the above kinematic range of 

a w->Xn = 1Q q ± 2 (stat.) ± 1.0(syst.) fib. 

The curves in Fig. ^|(a) are based on the effective flux of Eq. (g). For photoproduction, the x L 
dependence of the leading-neutron cross section in this model can be expressed as the integral 
over t (or p 2 T ) of Eq. (||): 



1 da 



LN 



a dx T 



oc (l-x.r^- 1 J f eS (x L ,t)dt, 



where it has been assumed that the jit total cross section has the same slow energy dependence 
as observed for other hadronic cross sections, characterised by the Pomeron intercept, «jp(0)) ~ 
1.08 |5T 
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The solid histogram in Fig. |8|(a) shows the result of a fit using the shape given by Eq. (Jjjfy plus a 
background term proportional to (1 —x L ). The fitQuses only statistical uncertainties and gives 
a good description of the data. Background terms of higher order in (1 — x L ) are not necessary. 

The HI collaboration also finds agreement between their leading-neutron data Q in the DIS 
regime and an OPE model. Simple leading-order OPE models (POMPYT and RAPGAP) also 
give a good description of the ZEUS data on dijets with a leading neutron in photoproduction 
|§. A recent NLO calculation |K| based on OPE is also in good agreement with the neutron- 
tagged dijet data. 



8.3 Vertex factorisation 



Leading baryon production in hadronic reactions has often been analysed in the context the 
Triple Regge expansion Hl6j|6l]| . If the pion trajectory dominates leading neutron production, 
as the data reported in this paper suggest, and if the si values are sufficiently high, then the 
7r-7r-Pomeron triple-Regge vertex will dominate. This is a reasonable approximation since s' is 
typically 80 GeV and factorisation should hold. The cross-section ratios are therefore given by 

d 2 a a ^ Xn (x L ,p 2 T ) I tfo^ Xn jx uV \) = aZ 
dx L dp 2 / dx L dp\ Otot' 

or, with a = 7 and b = p, 

1 d 2 a" /p ^ Xn (x L ,p 2 T ) 1 d 2 a pp ^ Xn (x L ,p 2 T ) 
c"tot dx L dp 2 T crfot dx L dp\ 



(13) 



To test this relation, the PHP data were integrated over p 2 T , assuming that the functional form 
of Eq. (|l|) also applies to photoproduction. The results for the PHP sample, expressed as the 
normalised cross section at p 2 T = 0, are shown in Fig. |8](b). These data are compared with ISR 
data 2£| from the reaction pp — > Xn. The shapes of the two distributions are similar but the 



PHP data are about a factor of two below the ISR pp results for x L > 0.3. These distributions, 
when integrated over all p\ and x L , represent the average number of neutrons per event, < n >. 
Using Eq. ( |T2"D to integrate over all p 2 T , the ZEUS PHP data yield < n >m 0.18 for the region 
x L > 0.28. The value 0.375±0.075 for 9 n < 150 mrad per hemisphere per inelastic pp collision 



was obtained at the ISR |27| and is a factor of two larger than the value obtained in the current 
analysis. 

The curve in Fig. |8](b), based on the OPE model with the effective flux of Eq. (|8]), gives a 
good representation of the shape of the ZEUS PHP data for x L > 0.6. In summary, while the 
shapes of the 'jp and pp distributions are similar, the cross-section scaling expressed in Eq. ( |T3"D 
is broken by about a factor of two. This breaking suggests that other Regge trajectories, with 
their associated interferences, are contributing differently to the photo- and hadroproduction 
reactions. 



5 The highest-x L bin was excluded from the fit because of the large systematic uncertainty arising from the 
energy scale. 
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8.4 Absorptive effects 



In 7*p scattering, the transverse size of the virtual photon decreases with increasing Q 2 , reducing 
the likelihood that the produced neutron rescatters on the hadronic component of the photon. 

At very high Q 2 , the virtual photon acts as a point-like probe. The HERA data thus offer 
the opportunity to study scattering with a target of fixed size (the proton) and a projectile of 
variable size (the photon). Figure |7](a) shows a clear reduction in the relative yield of neutrons in 
PHP and in the intermediate-Q 2 region compared to that in DIS at higher Q 2 . This reduction, 
which is displayed as a function of Q 2 in Fig. |9| for 0.64 < x L < 0.82, may be attributed to 
absorptive effects. Recently, a similar rise in the leading-proton yield as Q 2 is increased from 
the PHP region to Q 2 > 2.5 GeV 2 has been observed by the HI collaboration ||. 

The x L dependence of this effect is shown in Fig. [U], where the ratio -R abs (^L) is defined as 

r(Q 2 < 0.02 GeV 2 ,x L ) 



r(Q 2 > 4 GeV 2 , x L ) 



In order to minimise sensitivity to drifts in the relative energy scales, which particularly affects 
the higher values of x L , only PHP and DIS data taken simultaneously were used to evaluate the 



seven R ahE points highest in x L . Figure |X0| shows that R ahE is approximately constant, although 



some variation is observed for the higher x L values. 

The curves in Fig. [ll] are theoretical predictions based on OPE models p5| , |36| for the difference 
in absorption between Q 2 > 10 GeV 2 and photoproduction. The t dependence of the absorptive 
correction has been evaluated over a kinematic region similar to that of the present analysis to 
give the t-averaged correction. The calculations [^,|36| predict that absorption at the level of 
~10% may still be present for Q 2 > 10 GeV 2 . 

In conclusion, the characteristics of the production of leading neutrons in the range 0.64 < 
x L < 0.82, i.e. the rapid rise in rate with Q 2 and the saturation for the DIS data, are in broad 
agreement with the expectation from absorptive effects. 



8.5 Tests of limiting fragmentation and factorisation in x and Q 2 

The hypothesis of limiting fragmentation [|nj states that, in the high-energy limit, the cross 
section for the inclusive production of particle c in the fragmentation region of b in the reaction 
a + b — > c + X will approach a constant value. This suggestion, which was based on rather 



general geometrical arguments, was supported by measurements made at the ISR ||34|| . The 
expectation at HERA is that neutron production in the proton-fragmentation region of ep 
collisions will be independent of Q 2 and W |ll| . 



A more differential form of factorisation than that discussed in Section |8.3| states that the 
dependence of the cross section on the lepton variables (x and Q 2 ) should be independent of 
the baryon variables (x L and p 2 ). In terms of the structure function F2 N<4) , this can be written 

as 

F^ w (x 1 Q 2 ,x L} p T ) = f(x L ,p T ) ■ F(x,Q 2 ), 
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where / and F are arbitrary functions. The quantity _F 2 L N<3) is then given by 

F^ 3 \x,Q 2 ,x L ) = f~(x L )-F(x,Q 2 ), (14) 

where 

rp^(x L ) 

f(x L ) = / f(x L ,p T )dp T . 
Jo 

To investigate the W dependence, the data may be studied as a function of y, or x, for fixed Q 2 . 
Figure |TI| compares the cross-section ratio, r, as a function of Q 2 at fixed y by choosing three 
bins of y for a low (0.20 <x L < 0.64), a medium (0.64 <x L < 0.82) and a high (0.82 <x L < 1.0) 
x £ range. At low x L and fixed y, r rises with Q 2 , approximately doubling over the kinematic 
range of five orders of magnitude in Q 2 . As x L increases, the slope of r as a function of Q 2 
decreases and becomes negative in the high-x L range. 



Figure |T| shows the values of r for the intermediate-Q 2 data, separately for the low, medium 
and high ranges of x L . There is little variation of the ratio, r, as a function of y. The dotted 
horizontal lines in the figures, which show the average values of the points in each x L bin, are 
a satisfactory representation of the data. 

Figure shows r for the DIS data as a function of x in bins of Q 2 and x L . The same trends 



as in Fig. 11 are observed. At fixed Q 2 , there is a small, but systematic, dependence of r on x. 
At low x L , the ratio rises with x; for 0.64 < x L < 0.82, the ratio is approximately independent 
of x; at high x L , it falls with increasing x. The dotted lines show the results of fits of the form 
(xi/x) x , where X\ and A are fit parameters. Each x L range is fitted separately. 

These fits provide a good representation of the data except at the highest x L at high Q 2 . The 
data can also be fit to the form a + b In (x/xq) +cln {Q 2 /Qq), where a, b and c are fit parameters 
and Ql and xo are arbitrary scales. This gives a good representation of the data in all x L bins, 
as shown by the solid lines. 

In summary, the observed weak, logarithmic, variation as the size of the virtual photon changes 
represents a breakdown of limiting fragmentation and the factorisation as expressed in Eq. (|i~4"|). 



However, in the interval 0.64 < x L < 0.82, r for the DIS data is, to a good approximation, 
independent of the leptonic variables and is only a function of x L . For this x L range, Eq. ( p~0| ) 
implies that 

F 2 LN(3) (x,Q 2 ,x L )cxF 2 (a;,Q 2 ). (15) 



9 Determination of F 2 LN( ^ 

The neutron-tagged structure function, F 2 LN(3) , was determined using Eq. ([10|) . The F 2 (x,Q 2 
values were taken from the ZEUS Regge fit [p3] in the intermediate-Q 2 region, Q 2 < 1 GeV 2 
and from the ZEUS NLO QCD fit |3j for Q 2 > 1 GeV 2 . Figures |14| (intermediate-Q 2 data; 



and [T5| (DIS data) show _P 2 LN(3) for three representative ranges of x L as a function of x at fixed 
values of Q 2 . The normalisation of the three sets of fixed-x £ data varies because the range of 
p 2 T changes with x L when integrating over a fixed 9 n range. To guide the eye, the line shows 
the ZEUS Regge or NLO fit appropriately scaled by r in each x L bin divided by Ax L . The x 
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and Q 2 dependence of F 2 LN(3) is similar to that of F 2 , as expected since, as shown in Figs. |TT| 
and [13|, the ratio r has only a weak dependence on x and Q 2 at fixed x L . 

The HI Collaboration has published values of F% n in the kinematic region 2 < Q 2 < 50 
GeV 2 and 6 ■ 10~ 5 < x < 6 ■ 10~ 3 0. This region overlaps that covered by the ZEUS data 
for Q 2 > 4 GeV 2 . Although the HI data cover the ZEUS range in x L , the distribution of 
the neutrons is integrated only up to p™ ax = 0.2 GeV, so the ZEUS and HI data can only be 
directly compared at x L = 0.3, corresponding to p™ ax = 0.197 GeV. For higher x L , the ZEUS 
values must be reduced to account for the smaller p T range measured by HI. Figure |i6|(a) 
compares the ZEUS and HI values for F% ' at x L = 0.3 for the three bins of Q 2 where the 
data overlap. Figure |16|(b) shows the same comparison at x L = 0.7, where the ZEUS data have 
been adjusted to the transverse-momentum range p™ ax = 0.2 GeV using the form of Eq. fll2|) . 
The shapes of the distributions are in good agreement, as is the normalisation for the higher- 
x L selection. At low x L , although the HI values lie systematically above the corresponding 
ZEUS measurements, the HI uncertainties are highly correlated. In addition, contributions 
from photon showers, which populate the low x L region, have not been subtracted from the HI 
data. 



10 OPE and the pion structure function 

In a one-pion-exchange model, the structure function F 2 LN ^ can be written as 

F 2 LN(4) (x, Q 2 , x L , t) = U /p (x L7 t)F?(x/(l - x L ), Q\ t) (1 - A abs (Q 2 , x L , t)), (16) 

where j^i p is the flux of pions in the proton, F% is the structure function of the pion at a 
virtuality t, and A a b s is the absorptive correction before a t-integration is performed. 

There is now no simple factorisation of F 2 LN( ' 4 ' ) , although if F$ can be expressed as a power law 
in x-jr, where the Bjorken variable of the pion is defined as £„- = x/(l — x L ), as is the case for F 2 
of the proton with respect to x ||63|| , the approximate factorisation of Eq. ( |I~4D will be restored. 



Comparison of Eq. (14) with Eq. (|ToD (neglecting A) shows that if F is identified with F 2 , 
then r/Ax L may be identified with /. Given that r for the DIS data with 0.64 < x L < 0.82 
(where OPE dominates, see Section |10.1| ) is only a function of x L , as described in Section |8~5| , 



Eqs. (|T51) and (|IB|) imply that 

FZ( X7r ,Q 2 )<xF 2 (x,Q 2 ). 



10.1 Competing processes to OPE 



Several processes which compete with pion exchange as the mechanism for leading neutron 
production were ignored in Eq. (|T6|), namely: 



• diffractive dissociation in which the dissociated system decays to a state including a 
neutron 



16 



Diffractively produced events can be selected by requiring the presence of a large rapidity 
gap in the hadronic final state. For such events, the mass of the dissociated proton system 
is restricted to low values, M N < 4 GeV. 

An event is said to have a large rapidity gap (LRG) in the ZEUS detector if the pseudo- 
rapidity of the most-forward energy deposit with energy greater than 400 MeV (r? max ) is 
less than 1.8 |64]| . Figure |T7|(a) shows the 7? max distributions for both the neutron-tagged 
and inclusive DIS samples, where the latter has been normalised to the neutron-tagged 
sample for r] max > 1.8. For both r^ max < 1.8 and 7/ max > 1.8, the shape of the neutron-tagged 
distribution is similar to that of the inclusive distribution; however, there are relatively 
fewer LRG events in the neutron-tagged sample. The LRG events represent only 4% of 
the total number of DIS events with neutrons in the measured kinematic region, but are 
7% of the total number of DIS events. A reduction in the fraction of LRG events with 
a final-state neutron is expected since only proton diffractive dissociation or diffractive 
meson exchange (the Deck effect []S5|] ) can contribute. 

To investigate a possible x £ -dependence of the contribution of diffractive events, Fig. |17|(b) 
shows the ratio, Rlrg, of the neutron-tagged DIS events, selected by the LRG criterion, 
to all neutron-tagged DIS events, as a function of x L . The rise by a factor of three over 
the x L range shows that the LRG neutron-tagged events have a harder neutron energy 
spectrum than that of the inclusive neutron-tagged sample. It is clear that diffractive 
events are not a major source of leading neutrons at any value of x L . For the region 
0.64 <x L < 0.82, .Rlrg is 0.039 ± 0.001 (stat.); 

p and a 2 exchange 



Theoretical studies of neutron production in ep collisions ^ suggest that isovector 
exchanges other than the pion contribute less than 10% to neutron production at x L = 0.73 
and for the p T range of the present data. This is quite different than for leading proton 
production, where isoscalar Regge exchange provides the dominant contribution P,|T9||; 

isovector exchange leading to A production 

The p — > A(1236) transition, formed by 7r, p and a 2 exchange, can also contribute to 



neutron production fl~8|, ^3|,|66|-|68||. In this case, the neutron, which comes from the decay 



A — > nil or A + — > mr + , no longer has an energy determined by the energy of the 
exchanged meson. The neutron energy spectrum peaks near x L ~ 0.5 and extends only to 
x L ~ 0.7 |TJ|. It thus gives a small contribution in the 0.64 < x L < 0.82 bin. A comparison 



of the data on p — > n and p — > A ++ fl69|-[r2|j in charge-exchange reactions at Fermilab 



indicates that only about 6% of the forward neutrons come from the A channel. This 
observation agrees with theoretical estimates of the A7r contribution to the Fock state 



of the proton, which is approximately half that of nir []53|, p6|| . A calculation [67] shows 
that the contribution of p/a 2 exchange, plus the A contribution, to the hadronic charge- 
exchange reaction pp — > Xn could be as high as 30%. Since no analogous calculation 
exists for DIS, this only provides an indication of a possible background to the neutron 
production discussed in this paper; 

models other than one-particle exchange 

Monte Carlo studies, using standard DIS generators, show that these processes have 
a rate of neutron production a factor of three lower than the data and produce a neutron 
energy spectrum with the wrong shape, peaking at values of x L below 0.3. 
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In summary, the expectation for the processes listed above at (x L ) = 0.73 and {p T 2 ) = 0.08 GeV 2 
is that they contribute of the order of 20% of the leading neutron production. This estimate 
can be checked using the measured neutron-energy spectrum. The OPE fit to the differential 
cross-section da/dx L shown in Fig. ||(a) suggests that, at x L = 0.73, the residual background 
to OPE is < 10%, in reasonable accord with the studies discussed above. 



10.2 The pion structure function and the photon-pion total cross 
section 

Having established that the leading neutron data are dominated by pion exchange in the range 
0.64 < x L < 0.82, the OPE model can be used to determine the structure function of the pion, 
F£. The quantity F 2 LN ^ 3 ' 1 , defined in Section |2|, was obtained by integrating over t (or 

6 n ). The t dependence of F% is absorbed into the flux factor that describes the hadronic charge- 
exchange data [^4 -33]; the structure function of the real pion is then given by F£(x„, Q 2 ) = 



F^{x^Q\t = ml). 
Integrating Eq. (|T6D over t and rearranging leads to 

F*{x^Q 2 ) = T(Q 2 ,x L )F^ 3 \x,Q 2 ,x L ), (17) 

where T(Q 2 , x L ) is the inverse of the pion flux factor integrated over the measured t region and 
corrected for the t-averaged absorptive effect, <5 a b s : 

r(Q 2 ,x^ 1 



i - <Wx £ , Q 2 )) St~ U/ P (x L , t)dt 



As discussed in Section p\3| , the theoretical expectation |36| is that 5 a b s (x L , Q ) is less than 10% 
for Q 2 >10 GeV 2 . 

Thus, Eq. (|17D shows that, for the DIS region, F^x^.Q 2 ) is approximately proportional to 
F 2 LN (a?, Q 2 , x L ) for a fixed x L . In analogy to Eq. (|i~7l), the photon-pion total cross section can 
be written as 

aZ(wVT^T L ) = r(Q 2 = o, x L ) ■ d(jlp ~* Xn ^ Q 2 = °) (18) 

ax. 



L 



10.3 The pion-flux normalisation 

To determine the normalisation of the pion flux, and hence F% and a^t, a theoretical model is 
necessary for the form factor F(x L , t) in the effective flux factor of Eq. (0). Unfortunately, there 



is no consensus as to which of the competing models [|73|] is most appropriate. As examples 
of the extremes from the available range of the models, two possibilities are considered here. 
The evaluation of the normalisation is made at x L = 0.73, where the ratio r is approximately 
independent of x and Q 2 and OPE is expected to dominate. 

1. One possibility is to use the same flux as is employed in hadronic reactions. The data on 
the hadron charge-exchange reactions pn Xp, pp — > Xn and 7tn — > Xp p4]-p3| are in 
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rood agreement with OPE and the Bishari flux factor [JTJJ given by Eq. (Q). Using Eq. 
ffl), the cross section may be written as 



L- 



where a^t is the irp total cross section and the integration is over the kinematic range of 
the pn — > Xp measurement. The flux factor, / e g, can again be interpreted as an effective 
pion flux which already incorporates the effects of absorption, processes other than OPE 
and off-mass-shell effects. 

Assuming that all of these effects are the same for the photoproduction reaction 7p — * Xn 
as for the hadronic processes pn — > Xp and pp —>■ Xn then 

T(Q 2 = 0, (x L ) = 0.73) = ^^((^o.^^^ = ^75 = 2 ' 32 ' 
for the range 0.64 < x L < 0.82. 

Using Eq. (|T|) with da^ Xn /dx L evaluated from Fig. |(a) at x L = 0.73, a^ t (W = 
107 GeV) is 50 ± 2 fib. Assuming a Regge dependence such that a oc (W 2 ) apt <^~ 1 
with a F (0) = 1.08 |5TJ yields a^ t (W = 207 GeV) =56 ± 2 fib. Given that aZ(W = 
209 GeV) = 174 ± 13 fib ||, aZJa^ is then 0.32 ± 0.03. 



Using Fig. |9] to estimate the difference in the absorption for DIS and photoproduction as 
22% leads, for Q 2 > 4 GeV 2 , to 

FZ( X7r , Q 2 ) = 1.81 F 2 LN(3) (x, Q 2 , (x L ) = 0.73), 

or, combining Eqs. (|T0|) and (|17D, 



Q 2 ) = 10.05 • r(x, Q 2 , (x L ) = 0.73) ■ F 2 (x, Q 2 ). (19) 



2. As a second possibility, the additive quark model [[F4j] can be used, which predicts that, 
at x L = 0.73, 



-(1 — X L 
3 V 



,ajp(0)-l = q^q 



(20) 



<t(W) 

provided that cr^t and a^ t have the same energy dependence at high energies, governed 
by the Pomeron intercept, ajp(0) = 1.08. Combining Eq. ( |TT] ) with Eqs. (|TJ) and fl2Tf ) 
and using Table |3| gives 

r(Q> = 0, W = 0.73) = |(1 - 0.73)^-. -_^__ = JL52_ = 4 .8. 

Using i? aba to adjust for the difference in absorptive effects between photoproduction and 
DIS, thereby reducing T, yields, for Q 2 > 4 GeV 2 : 



F^x^Q 2 ) = 3.74F 2 LN(3) (x,Q 2 , (x L ) = 0.73) 
= 20.78 • r(x, Q 2 , (x L ) = 0.73) • F 2 (x, Q 2 ). 



(21) 
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In summary, methods 1 and 2 differ in normalisation by about a factor of two. This is ap- 
proximately the same factor that was observed in Fig. |8|(b) when comparing the neutron rates 
presented here and the pp data from the ISR. 

Other methods give factors between those described above. For example, models of the pion 
flux factor based on d/u data give form factors, F(x L , t) < 1 |73| , [75| and, therefore, a value of T 
that is higher than that given by method 1. Normalisation of F£ using existing measurements 
at high x w |38|-f44|l, as parameterised, for example, by Gliick, Reya and Vogt 
value that is lower than that from method 2. 



T6fl , gives a T 



10.4 Results and comparison with models 



Figure [L8| shows F 2 n,EF (x 7T ,Q 2 ), which is F^x^^Q 2 ) as a function of x^ evaluated from Eq. 
(|19D. This evaluation has been made for the range 0.64 < x L < 0.82 where OPE is expected 
to dominate and where the fraction of events with a leading neutron is observed to be approx- 
imately independent of x and Q 2 . Evaluations for nearby ranges of x L give consistent results. 
Although the normalisation of F% is uncertain to a factor of two, the x and Q 2 dependence 
is well measured. The dotted lines show the shape of F 2 of the proton, scaled by 0.361. It is 
striking that F£ has approximately the same x and Q 2 dependence as F 2 of the proton. 



Figure ECU] shows F^ A( ^ M using the additive quark model to normalise the cross section. The 



comparison of these results with the solid curves, which represent the absolute predictions of 
Gliick, Reya and Vogt (GRV) [76), shows that the normalisation is such that the curves are 
closer to the data than when using the effective flux shown in Fig. [18], although the data are 
still higher than the predictions, particularly at the lower values of Q 2 . 

There have been several recent attempts to understand the quark structure of the pion on the 
basis of simple theoretical ideas: 

1. In the model of Gliick, Reya and Vogt [ 76H , the pion sea is generated dynamically by the 



QCD evolution equations from the valence-quark and valence-gluon distributions of the 
pion at a low scale of 0.3 GeV 2 . The valence gluon is, by assumption, proportional to the 
valence quark. The valence quarks are taken from fits to experimental data [|3"8"|-|4"4|| at 
high x n . This model qualitatively describes the shape of the measurement of F£ as seen 
in Figs ITS and O; 



The prediction for F£ from Sutton et al. (SMRS) [77 is also shown in Figure |T^. This 



analysis is an NLO QCD fit to Drell-Yan and prompt-photon data from fixed-target 
experiments. 

Since the data corresponds to x n > 0.2, the predictions for the sea region are uncertain. 
However, the SMRS fits, when used in an OPE Monte Carlo simulation, gave a good 
description of the ZEUS data on leading neutrons plus dijets in photoproduction [[|. The 
latter correspond to the region x n > 10~ 2 , whereas the present data extend to much lower 
x n and show a very different shape than the F£ from the SMRS parameterisation; 



3. In the constituent model of Altarelli et al J75|] , the structure function of any hadron 
is determined by that of its constituent quarks, and, in particular, the pion structure 



20 



function can be predicted from the known nucleon structure function [|79|| . The structure 
function of a hadron h, F£, at low x, where the contribution of the constituent valence 
quarks is negligible, is given by 

F?(x,Q 2 )= f 1 f h (z)F«(x/z,Q 2 )dz, 

J x 

where fh{ z ) is the density of quarks carrying momentum fraction z of the hadron h, and 
F| is the structure function of the valence quarks. It is assumed that, at low x, the 
structure function F$ is independent of the flavour. If, as x — > 0, F$(x) ~ x~~ x ^ \ where 
X(Q 2 ) is independent of x [|63[| , all of the F^ structure functions are proportional to the 
same function at low x. This prediction is in agreement with the present data since 
F£ oc Fi- This model gives a normalisation for F£ close to p^ A Q M -^ 

4. Nikolaev, Speth and Zoller have explicitly studied the small-x behavior of F% using a 
colour-dipole BFKL-Regge expansion ||80|| . They find that a good approximation is: 



F-(*,Q 2 )^Ff (|r,Q 2 ) 



Since F2 has a power-law behaviour in x for the ZEUS kinematic region, the predicted 
shape is in reasonable agreement with the data. 

Thus, the data presented in this paper provide additional constraints on the shape of the pion 
structure function for values of x^ < 10 -2 . 



11 Conclusions 

Leading neutron production for x L > 0.2 and 9 n < 0.8 mrad has been studied in neutral 
current ep collisions at HERA from photoproduction to Q 2 ~ 10 4 GeV 2 . The cross section for 
the production of leading neutrons has been determined as a ratio r relative to the inclusive 
neutral current cross section, thereby reducing the systematic uncertainty considerably. 

For the above angular range, which defines a transverse-momentum acceptance of p T < 0.656 x L 
GeV, the neutron energy spectrum exhibits a broad peak at x L ~ 0.75. It is approximately one 
third of its maximum at x L = 0.2, and falls to zero approaching the kinematic limit at x L = 1. 
The rate of neutron production in the measured phase space in photoproduction is about half 
the rate observed in hadroproduction and thus the simplest form of vertex factorisation is 
broken. A comparison of the neutron yield in the three kinematic regions of photoproduction, 
intermediate Q 2 and DIS shows an increase of about 20% between photoproduction and DIS, 
saturating for Q 2 > 4 GeV 2 . This can be attributed to the decrease in absorptive effects as the 
transverse size of the incident photon decreases. 

Generally, there is no strong dependence of the ratio r on x and Q 2 ; however, at low x L and 
fixed y, the ratio r rises with Q 2 , approximately doubling over the kinematic range of five orders 
of magnitude in Q 2 . For the high-x £ range, there is a tendency for r to decrease at fixed y. 
Thus, limiting fragmentation and factorisation do not hold. 
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The structure-function F 2 has been measured over three orders of magnitude in x and Q 2 
and for the range 0.2 < x L < 0.8 with 8 n < 0.8 mrad. The structure function of the pion, 
F£, has been extracted, up to uncertainties in the overall normalisation, in the framework of a 
one-pion-exchange model in the range 0.64 < x L < 0.82, where one-pion-exchange dominates. 
It has approximately the same x and Q 2 dependence as F 2 of the proton. The data provide 
new constraints on the shape of the pion structure function for x n < 10~ 2 . 
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Appendix I: Tables of relative cross sections 

This appendix presents tables of the relative cross sections for leading neutron production in 
neutral current interactions, ep — > e'Xn. The cross sections are measured relative to the total 
inclusive neutral current cross section, ep — > e'X. 

The x L bins, the acceptance and the acceptance and energy-scale uncertainties for each bin are 
given in Table III. Only for 0.88 < x L < 1 is the centroid of the bin, 0.92, significantly different 
from the bin centre, 0.94. 

The values of b, the exponential slope of the p\ distribution, are given in Table |2|. 

Table [| gives the relative differential cross section in x L , {l/a){da hN /dx L ), for leading neutron 
production integrated over x and Q 2 for the three kinematic regimes: a) photoproduction 
(Q 2 < 0.02 GeV 2 ), (b) the intermediate-Q 2 region (0.1 < Q 2 < 0.74 GeV 2 ), and (c) DIS 
(Q 2 > 4 GeV 2 ). 

A change in the energy scale represents a stretch or compression (dilation) of the energy distri- 
bution. Although the dilation is only ±2%, it becomes the dominant normalisation systematic 
uncertainty at high x L , where the energy spectrum is falling quickly (Table |^). 

The next twelve tables give the fraction of events with a leading neutron in DIS (Q 2 > 4 GeV 2 ) 
in bins of Q 2 and x, while the final three tables give the fraction of events with a leading neutron 
in the intermediate-Q 2 region (0.1 < Q 2 < 0.74 GeV 2 ) in bins of Q 2 and y. 
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range 


<x L > 


Acceptance 
(%) 


Acceptance 
uncertainty 

(%) 


Energy-scale 
uncertainty contribution 

(%) 


0.20 - 0.28 


0.24 


21.5 


+2/-2 


-5/+4 


0.28 - 0.34 


0.31 


21.2 


+2/-2 


-5/+4 


0.34 - 0.40 


0.37 


21.7 


+2/-3 


-5/+4 


0.40 - 0.46 


0.43 


22.2 


+3/-3 


-5/+4 


0.46 - 0.52 


0.49 


22.5 


+3/-3 


-4/+5 


0.52 - 0.58 


0.55 


21.9 


+4/-3 


-4/+5 


0.58 - 0.64 


0.61 


23.0 


+4/-4 


-4/+4 


0.64 - 0.70 


0.67 


22.9 


+5/-5 


-4/+4 


0.70 - 0.76 


0.73 


23.3 


+5/-5 


-2/+2 


0.76 - 0.82 


0.79 


23.7 


+5/-6 


+2/-1 


0.82 - 0.88 


0.85 


25.4 


+6/-6 


+5/-8 


0.88 - 1.00 


0.92 


32.1 


+6/-6 


+25/-26 



Table 1: The x L bins, their acceptance, and the acceptance uncertainty. The right hand 
column shows the contribution from the energy-scale uncertainty. Note that the energy-scale 
uncertainty is completely correlated between bins. 
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b 

GeV- 2 


±(stat.) 
GeV- 2 


0.49 


4.10 


1.10 


0.54 


4.47 


1.00 


0.58 


3.80 


1.10 


0.62 


6.33 


1.00 


0.65 


8.12 


0.90 


0.69 


8.56 


0.90 


0.73 


6.06 


0.80 


0.76 


9.46 


0.90 


0.80 


8.78 


0.80 


0.84 


12.81 


1.10 


0.88 


9.65 


0.90 


0.95 


8.63 


1.10 



Table 2: The values of the b slopes and their statistical errors in bins of x L for the 1995 DIS 
data. 
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X L 

range 


<x L > 


Q 2 < 0.02 GeV 2 
meas. ± stat. 


0.1 < Q 2 < 0.74 GeV 2 
meas. ± stat. 

(%) 


Q 2 > 4 GeV 2 
meas. ± stat. 

(%) 


0.20 - 0.28 


0.24 


4.12 ± 0.60 


4.70 ± 0.36 


5.57 ± 0.13 


0.28 - 0.34 


0.31 


4.40 ± 0.71 


6.37 ± 0.46 


7.19 ± 0.17 


0.34 - 0.40 


0.37 


4.64 ± 0.72 


5.83 ± 0.42 


8.51 ± 0.18 


0.40 - 0.46 


0.43 


6.31 ± 0.84 


7.66 ± 0.50 


9.44 ± 0.19 


0.46 - 0.52 


0.49 


6.90 ± 0.87 


7.72 ± 0.49 


10.38 ± 0.20 


0.52 - 0.58 


0.55 


7.94 ± 0.34 


10.64 ± 0.57 


12.38 ± 0.22 


0.58 - 0.64 


0.61 


9.40 ± 0.34 


11.03 ± 0.59 


13.43 ± 0.22 


0.64 - 0.70 


0.67 


11.15 ± 0.38 


12.69 ± 0.62 


15.28 ± 0.24 


0.70 - 0.76 


0.73 


12.51 ± 0.41 


14.06 ± 0.64 


15.83 ± 0.24 


0.76 - 0.82 


0.79 


12.15 ± 0.39 


14.41 ± 0.67 


14.85 ± 0.23 


0.82 - 0.88 


0.85 


9.18 ± 0.31 


11.11 ± 0.56 


11.66 ± 0.20 


0.88 - 1.00 


0.92 


2.59 ± 0.09 


3.25 ± 0.19 


3.74 ± 0.07 



Table 3: Values in percent of (1 / a) (da LN / dx L ) = dr/dx L for leading neutrons with 9 n < 0.8 
mrad. For the PHP data, Q 2 < 0.02 GeV 2 , the normalisation uncertainty is ±5%, and for 
x L > 0.52 there is a further normalisation uncertainty of ±4%. The statistical error for the 
PHP data with x L > 0.52 includes a contribution from the uncertainty in the FNC-trigger 
efficiency For the intermediate-Q 2 and DIS data, the normalisation uncertainty is ±4%. 
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0.2 <x L < 0.28 


Q 2 

GeV 2 


Q 2 

range 


%BJ 


%BJ 

range 


ratio (%) 
meas. ± stat 


7 


4 - 10 


1.1 -10^ 4 

2.1 -10" 4 

4.2 -10" 4 
8.5 -10" 4 
17.0 -10" 4 


8.0- 15.0 -lO" 5 
1.5 - 3.0 -10" 4 
3.0 - 6.0 -10" 4 

6.0 - 12.0 -10" 4 
1.2 - 2.4 -10- 3 


0.48 ± 0.04 
0.40 ± 0.03 
0.43 ± 0.03 
0.46 ± 0.04 
0.50 ± 0.05 


15 


10 - 20 


2.1 -10- 4 

4.2 -10" 4 
8.5 -10" 4 
17.0 -10" 4 
49.0 -10" 4 


1.5 - 3.0 -10^ 4 
3.0 - 6.0 -10" 4 

6.0 - 12.0 -10" 4 
1.2 - 2.4 -10~ 3 

2.4 - 10.0 -lO" 3 


0.40 ± 0.05 
0.39 ± 0.03 
0.42 ± 0.03 
0.43 ± 0.04 
0.42 ± 0.03 


30 


20 - 40 


4.2 -lO" 4 
8.5 -10^ 4 
17.0 -HT 4 
49.0 -lO" 4 


3.0 - 6.0 -10^ 4 
6.0 - 12.0 -10- 4 

1.2 - 2.4 -10- 3 
2.4 - 10.0 -10~ 3 


0.55 ± 0.06 
0.50 ± 0.04 
0.43 ± 0.03 
0.49 ± 0.03 


60 


40 - 80 


8.5 -10~ 4 
17.0 -10" 4 
49.0 -10" 4 
32.0 -HT 3 


6.0 - 12.0 -10- 4 
1.2 - 2.4 -10- 3 
2.4 - 10.0 -lO^ 3 
1.0 - 10.0 -lo- 2 


0.49 ± 0.08 
0.46 ± 0.05 
0.39 ± 0.03 
0.42 ± 0.04 


120 


80 - 160 


17.0 -lO" 4 
49.0 -10" 4 
32.0 -10- 3 


1.2 - 2.4 -10- 3 
2.4 - 10.0 -lO" 3 
1.0 - 10.0 -lO" 2 


0.53 ± 0.12 
0.42 ± 0.05 
0.53 ± 0.05 


240 


160 - 320 


49.0 -10" 4 
32.0 -10- 3 


2.4- 10.0 -10- 3 
1.0 - 10.0 -lo- 2 


0.42 ± 0.09 
0.56 ± 0.07 


480 


320 - 640 


32.0 -lO" 3 


1.0 - 10.0 -lo- 2 


0.44 ± 0.09 


1000 


640 - 10000 


32.0 -lO" 3 


1.0 - 10.0 -10^ 2 


0.37 ± 0.11 



Acceptance uncertainty: +2/-2 % 



Energy scale uncertainty of ±2 %: -5/+4 % 
Normalization error: 4 % 



Table 4: The fraction of events with a leading neutron at x L = 0.24 in bins of x and Q 2 . 
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0.28 <x L < 0.34 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


io- 4 


8.0 - 


15.0 


•10 


-5 


0.34 ± 0.04 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.37 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.41 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.50 ± 0.04 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.46 ± 0.04 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.40 ± 0.05 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.40 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.45 ± 0.03 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.38 ± 0.04 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.45 ± 0.04 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.47 ± 0.06 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.49 ± 0.04 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.51 ± 0.04 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.50 ± 0.03 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.37 ± 0.07 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.38 ± 0.04 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.52 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.51 ± 0.05 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


10" 


3 


0.24 ± 0.08 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.43 ± 0.05 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.56 ± 0.05 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.43 ± 0.09 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.40 ± 0.06 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.55 ± 0.11 


1000 


640 - 


L0000 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.45 ± 0.12 



Acceptance uncertainty: +2/-2 % 



Energy scale uncertainty of ±2 %: -5/+4 % 
Normalization error: 4 % 



Table 5: The fraction of events with a leading neutron at xl = 0.31 in bins of x and Q 2 . 
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0.34 <x L < 0.4 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


io- 4 


8.0 - 


15.0 


•10 


-5 


0.45 ± 0.04 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.46 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.52 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.50 ± 0.04 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.55 ± 0.05 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.42 ± 0.05 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.52 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.54 ± 0.04 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.48 ± 0.04 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.59 ± 0.04 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.45 ± 0.06 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.49 ± 0.04 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.55 ± 0.04 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.58 ± 0.03 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.56 ± 0.09 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.53 ± 0.05 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.51 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.58 ± 0.05 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


io- 


3 


0.52 ± 0.12 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.51 ± 0.05 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.54 ± 0.05 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.52 ± 0.10 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.59 ± 0.07 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.63 ± 0.11 


1000 


640 - 


L0000 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.54 ± 0.13 



Acceptance uncertainty: +2/-3 % 



Energy scale uncertainty of ±2 %: -5/+4 % 
Normalization error: 4 % 



Table 6: The fraction of events with a leading neutron at xl = 0.37 in bins of x and Q 2 . 
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0.4 <x L < 0.46 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


10~ 4 


8.0 - 


15.0 


•10 


-5 


0.51 ± 0.04 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.55 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.55 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.55 ± 0.04 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.62 ± 0.05 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.52 ± 0.05 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.55 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.55 ± 0.04 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.66 ± 0.05 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.60 ± 0.04 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.41 ± 0.05 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.46 ± 0.03 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.54 ± 0.04 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.65 ± 0.03 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.48 ± 0.08 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.64 ± 0.06 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.53 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.67 ± 0.06 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


10" 


3 


0.65 ± 0.14 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.59 ± 0.06 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.64 ± 0.06 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.77 ± 0.12 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.70 ± 0.08 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.77 ± 0.12 


1000 


640 - 


L0000 


32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.75 ± 0.16 



Acceptance uncertainty: +3/-3 % 



Energy scale uncertainty of ±2 %: -5/+4 % 
Normalization error: 4 % 



Table 7: The fraction of events with a leading neutron at xl = 0.43 in bins of x and Q 2 . 
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0.46 <x L < 0.52 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


10~ 4 


8.0 - 


15.0 


•10 


-5 


0.51 ± 0.04 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.59 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.57 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.67 ± 0.05 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.69 ± 0.05 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.55 ± 0.05 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.59 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.65 ± 0.04 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.67 ± 0.05 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.62 ± 0.04 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.50 ± 0.06 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.51 ± 0.04 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.66 ± 0.04 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.71 ± 0.03 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.55 ± 0.09 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.68 ± 0.06 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.69 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.80 ± 0.06 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


10" 


3 


0.92 ± 0.16 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.64 ± 0.06 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.86 ± 0.07 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.58 ± 0.10 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.69 ± 0.08 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.74 ± 0.12 


1000 


640 - 


L0000 


32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.81 ± 0.16 



Acceptance uncertainty: +3/-3 % 



Energy scale uncertainty of ±2 %: -4/+5 % 
Normalization error: 4 % 



Table 8: The fraction of events with a leading neutron at x L = 0.49 in bins of x and Q 2 . 
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0.52 < x L < 0.58 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


io- 4 


8.0 - 


15.0 


•10 


-5 


0.71 ± 0.05 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.70 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.66 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.73 ± 0.05 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.76 ± 0.06 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.71 ± 0.06 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.76 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.79 ± 0.04 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.72 ± 0.05 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.84 ± 0.05 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.72 ± 0.07 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.72 ± 0.04 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.79 ± 0.05 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.81 ± 0.04 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.75 ± 0.10 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.69 ± 0.06 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.71 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.88 ± 0.06 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


10" 


3 


0.77 ± 0.15 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.80 ± 0.07 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.88 ± 0.07 


240 


160 - 


320 


49.0 


•IO" 4 


2.4- 


10.0 


•10 


-3 


0.83 ± 0.12 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.97 ± 0.09 


480 


320 - 


640 


32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.88 ± 0.13 


1000 


640 - 


L0000 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.86 ± 0.17 



Acceptance uncertainty: +4/-3 % 



Energy scale uncertainty of ±2 %: -4/+5 % 
Normalization error: 4 % 



Table 9: The fraction of events with a leading neutron at xl = 0.55 in bins of x and Q 2 . 
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0.58 <x L < 0.64 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


io- 4 


8.0 - 


15.0 


•10 


-5 


0.78 ± 0.05 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.79 ± 0.04 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.81 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.78 ± 0.05 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.88 ± 0.06 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.68 ± 0.06 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.77 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.78 ± 0.04 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.85 ± 0.05 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.86 ± 0.05 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.76 ± 0.07 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.80 ± 0.04 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.75 ± 0.04 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.84 ± 0.04 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.69 ± 0.10 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.83 ± 0.06 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.87 ± 0.05 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.87 ± 0.06 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


10" 


3 


0.52 ± 0.12 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.80 ± 0.06 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.83 ± 0.06 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.84 ± 0.12 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


1.04 ± 0.09 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.89 ± 0.13 


1000 


640 - 


L0000 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.82 ± 0.16 



Acceptance uncertainty: +4/-4 % 



Energy scale uncertainty of ±2 %: -4/+4 % 
Normalization error: 4 % 



Table 10: The fraction of events with a leading neutron at x L = 0.61 in bins of x and Q 2 . 
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0.64 <x L < 0.7 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


io- 4 


8.0 - 


15.0 


•10 


-5 


0.86 ± 0.06 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.91 ± 0.04 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.92 ± 0.05 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


1.00 ± 0.06 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.85 ± 0.06 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.93 ± 0.07 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.93 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.91 ± 0.05 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.97 ± 0.05 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.92 ± 0.05 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.97 ± 0.08 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.86 ± 0.05 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.86 ± 0.05 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.95 ± 0.04 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.89 ± 0.11 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.81 ± 0.06 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.92 ± 0.05 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.97 ± 0.07 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


io- 


3 


0.82 ± 0.15 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.91 ± 0.07 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.95 ± 0.07 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


1.04 ± 0.13 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.90 ± 0.09 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


1.10 ± 0.14 


1000 


640 - 


L0000 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.98 ± 0.18 



Acceptance uncertainty: +5/-5 % 



Energy scale uncertainty of ±2 %: -4/+4 % 
Normalization error: 4 % 



Table 11: The fraction of events with a leading neutron at x L = 0.67 in bins of x and Q 2 . 
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0.7 <x L < 0.76 


Q 2 

GeV 2 


Q 2 

range 


%BJ 


%BJ 

range 


ratio (%) 
meas. ± stat 


7 


4 - 10 


1.1 -10^ 4 

2.1 -10" 4 

4.2 -10" 4 
8.5 -10" 4 
17.0 -10" 4 


8.0- 15.0 -lO" 5 
1.5 - 3.0 -10" 4 
3.0 - 6.0 -10" 4 

6.0 - 12.0 -10" 4 
1.2 - 2.4 -10- 3 


0.95 ± 0.06 
0.90 ± 0.04 
0.96 ± 0.05 
0.98 ± 0.05 
0.92 ± 0.06 


15 


10 - 20 


2.1 -10- 4 

4.2 -10" 4 
8.5 -10" 4 
17.0 -10" 4 
49.0 -10" 4 


1.5 - 3.0 -10^ 4 
3.0 - 6.0 -10" 4 

6.0 - 12.0 -10" 4 
1.2 - 2.4 -10~ 3 

2.4 - 10.0 -lO" 3 


0.84 ± 0.07 
0.88 ± 0.04 
0.94 ± 0.05 
1.02 ± 0.06 
1.09 ± 0.05 


30 


20 - 40 


4.2 -lO" 4 
8.5 -10^ 4 
17.0 -HT 4 
49.0 -lO" 4 


3.0 - 6.0 -10^ 4 
6.0 - 12.0 -10- 4 

1.2 - 2.4 -10- 3 
2.4 - 10.0 -10~ 3 


0.85 ± 0.08 
0.83 ± 0.04 
0.93 ± 0.05 
0.96 ± 0.04 


60 


40 - 80 


8.5 -10~ 4 
17.0 -10" 4 
49.0 -10" 4 
32.0 -HT 3 


6.0 - 12.0 -10- 4 
1.2 - 2.4 -10- 3 
2.4 - 10.0 -10- 3 
1.0 - 10.0 -lo- 2 


1.03 ± 0.12 
0.87 ± 0.06 
0.91 ± 0.05 
1.06 ± 0.07 


120 


80 - 160 


17.0 -lO" 4 
49.0 -10" 4 
32.0 -10- 3 


1.2 - 2.4 -10- 3 
2.4 - 10.0 -lO" 3 
1.0 - 10.0 -10~ 2 


1.11 ± 0.17 
0.83 ± 0.06 
0.94 ± 0.07 


240 


160 - 320 


49.0 -10" 4 
32.0 -10- 3 


2.4- 10.0 -10- 3 
1.0 - 10.0 -lo- 2 


0.96 ± 0.13 
1.01 ± 0.09 


480 


320 - 640 


32.0 -lO" 3 


1.0 - 10.0 -lo- 2 


1.10 ± 0.14 


1000 


640 - 10000 


32.0 -lO" 3 


1.0 - 10.0 -10^ 2 


1.16 ± 0.19 



Acceptance uncertainty: +5/-5 % 



Energy scale uncertainty of ±2 %: -2/+2 % 
Normalization error: 4 % 



Table 12: The fraction of events with a leading neutron at x L = 0.73 in bins of x and Q 2 . 
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0.76 <x L < 0.82 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


io- 4 


8.0 - 


15.0 


•10 


-5 


0.98 ± 0.06 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.83 ± 0.04 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.90 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.86 ± 0.05 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.93 ± 0.06 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.94 ± 0.07 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.89 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.90 ± 0.05 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.83 ± 0.05 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.87 ± 0.05 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.77 ± 0.07 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.89 ± 0.05 








17.0 


•io~ 4 


1.2 - 


2.4 


10" 


3 


0.82 ± 0.05 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.89 ± 0.04 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.91 ± 0.11 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.83 ± 0.06 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.89 ± 0.05 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.98 ± 0.07 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


10" 


3 


1.22 ± 0.18 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.67 ± 0.06 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.89 ± 0.07 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.93 ± 0.12 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.84 ± 0.08 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.60 ± 0.10 


1000 


640 - 


L0000 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.89 ± 0.17 



Acceptance uncertainty: +5/-6 % 



Energy scale uncertainty of ±2 %: +2/-1 % 
Normalization error: 4 % 



Table 13: The fraction of events with a leading neutron at x L = 0.79 in bins of x and Q 2 . 
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0.82 <x L < 0.88 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


10~ 4 


8.0 - 


15.0 


•10 


-5 


0.73 ± 0.05 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.71 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.65 ± 0.04 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.74 ± 0.05 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.74 ± 0.05 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.72 ± 0.06 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.70 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.70 ± 0.04 








17.0 


•HT 4 


1.2 - 


2.4 


10" 


3 


0.72 ± 0.04 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.75 ± 0.04 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.73 ± 0.07 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.65 ± 0.04 








17.0 


•io- 4 


1.2 - 


2.4 


10" 


3 


0.64 ± 0.04 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.68 ± 0.03 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.61 ± 0.09 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.69 ± 0.05 








49.0 


•10~ 4 


2.4- 


10.0 


•10 


-3 


0.61 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.66 ± 0.05 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


io- 


3 


0.72 ± 0.13 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.63 ± 0.05 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.60 ± 0.05 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.43 ± 0.08 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.40 ± 0.05 


480 


320 - 


640 


32.0 


•10~ 3 


1.0 - 


10.0 


•10 


-2 


0.40 ± 0.08 


1000 


640 - 


L0000 


32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.57 ± 0.13 



Acceptance uncertainty: +6/-6 % 



Energy scale uncertainty of ±2 %: +5/-8 % 
Normalization error: 4 % 



Table 14: The fraction of events with a leading neutron at x L = 0.85 in bins of x and Q 2 . 
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0.88 < x L < 1 



Q 2 


Q 


2 


X 


B.I 










ratio (%) 


GeV 2 


ran 


ge 

o 








range 




meas. ± stat 


7 


4 - 


10 


1.1 


10~ 4 


8.0 - 


15.0 


•10 


-5 


0.52 ± 0.04 








2.1 


10" 4 


1.5 - 


3.0 


10" 


4 


0.53 ± 0.03 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.48 ± 0.03 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.41 ± 0.03 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.45 ± 0.04 


15 


10 - 


20 


2.1 -IO" 4 


1.5 - 


3.0 


10" 


4 


0.51 ± 0.04 








4.2 


10" 4 


3.0 - 


6.0 


10" 


4 


0.45 ± 0.02 








8.5 


10" 4 


6.0 - 


12.0 


•10 


-4 


0.44 ± 0.03 








17.0 


•HT 4 


1.2 - 


2.4 


io- 


3 


0.43 ± 0.03 








49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.38 ± 0.03 


30 


20 - 


40 


4.2 


10~ 4 


3.0- 


6.0 


10" 


4 


0.40 ± 0.04 








8.5 


10~ 4 


6.0 - 


12.0 


•10 


-4 


0.39 ± 0.03 








17.0 


•io- 4 


1.2 - 


2.4 


10" 


3 


0.41 ± 0.03 








49.0 


•10" 4 


2.4 - 


10.0 


•10 


-3 


0.43 ± 0.02 


60 


40 - 


80 


8.5 


10~ 4 


6.0- 


12.0 


•10 


-4 


0.52 ± 0.07 








17.0 


•10~ 4 


1.2 - 


2.4 


10" 


3 


0.34 ± 0.03 








49.0 


•IO" 4 


2.4- 


10.0 


•10 


-3 


0.38 ± 0.03 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.34 ± 0.03 


120 


80 - 


160 


17.0 


•10" 4 


1.2 - 


2.4 


io- 


3 


0.45 ± 0.09 








49.0 


•io- 4 


2.4 - 


10.0 


•10 


-3 


0.36 ± 0.04 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.35 ± 0.04 


240 


160 - 


320 


49.0 


•10" 4 


2.4- 


10.0 


•10 


-3 


0.41 ± 0.07 








32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.22 ± 0.04 


480 


320 - 


640 


32.0 


•IO" 3 


1.0 - 


10.0 


•10 


-2 


0.25 ± 0.06 


1000 


640 - 


L0000 


32.0 


■10~ 3 


1.0 - 


10.0 


•10 


-2 


0.20 ± 0.07 



Acceptance uncertainty: +6/-6 % 



Energy scale uncertainty of ±2 %: +25/-26 % 
Normalization error: 4 % 



Table 15: The fraction of events with a leading neutron at x L = 0.94 in bins of x and Q 2 . 
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0.2 < x L < 0.64 



Q 2 


Q 2 


y 


y 




ratio (%) 


GeV 2 


range 




ran 


ge 


meas. ± stat 


0.11 


0.10 - 0.13 


0.60 


0.54- 


0.64 


3.15 ± 0.43 






0.70 


0.64- 


0.74 


2.85 ± 0.39 


0.15 


0.13 - 0.17 


0.40 


0.37- 


0.45 


3.72 ± 0.47 






0.50 


0.45 - 


0.54 


2.76 ± 0.34 






0.60 


0.54- 


0.64 


3.49 ± 0.39 






0.70 


0.64 - 


0.74 


2.76 ± 0.41 


0.20 


0.17 - 0.21 


0.26 


0.23 - 


0.30 


3.23 ± 0.41 






0.33 


0.30 - 


0.37 


3.09 ± 0.38 






0.40 


0.37- 


0.45 


3.41 ± 0.38 






0.50 


0.45 - 


0.54 


3.70 ± 0.43 






0.60 


0.54 - 


0.64 


2.96 ± 0.43 






0.70 


0.64- 


0.74 


3.98 ± 0.63 


0.25 


0.21 - 0.27 


0.20 


0.16 - 


0.23 


3.23 ± 0.28 






0.26 


0.23 - 


0.30 


3.46 ± 0.31 






0.33 


0.30 - 


0.37 


3.24 ± 0.33 






0.40 


0.37- 


0.45 


4.04 ± 0.39 






0.50 


0.45 - 


0.54 


3.28 ± 0.39 






0.60 


0.54- 


0.64 


3.83 ± 0.49 


0.30 


0.27 - 0.35 


0.12 


0.08 - 


0.16 


3.88 ± 0.29 






0.20 


0.16 - 


0.23 


3.82 ± 0.30 






0.26 


0.23 - 


0.30 


3.03 ± 0.30 






0.33 


0.30 - 


0.37 


2.86 ± 0.32 






0.40 


0.37- 


0.45 


3.36 ± 0.40 






0.50 


0.45 - 


0.54 


2.97 ± 0.41 


0.40 


0.35 - 0.45 


0.12 


0.08 - 


0.16 


3.05 ± 0.28 






0.20 


0.16 - 


0.23 


3.02 ± 0.30 






0.26 


0.23 - 


0.30 


3.61 ± 0.38 






0.33 


0.30 - 


0.37 


3.94 ± 0.45 






0.40 


0.37- 


0.45 


4.45 ± 0.52 


0.50 


0.45 - 0.58 


0.12 


0.08 - 


0.16 


3.35 ± 0.33 






0.20 


0.16 - 


0.23 


3.02 ± 0.34 






0.26 


0.23 - 


0.30 


2.50 ± 0.36 


0.65 


0.58 - 0.74 


0.20 


0.16 - 


0.23 


4.01 ± 0.54 






0.26 


0.23 - 


0.30 


3.46 ± 0.68 



Acceptance uncertainty: +2/-2 % 



Energy scale uncertainty of ±2 %: -4/+4 % 
Normalization error: 4 % 



Table 16: The fraction of events with a leading neutron at Xl = 0.42 in bins of y and Q 2 for 
the BPC region. 
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0.64 < x L < 0.82 



Q 2 


Q 2 


y 


y 




ratio (%) 


GeV 2 


range 




ran 


ge 


meas. ± stat 


0.11 


0.10 - 0.13 


0.60 


0.54- 


0.64 


2.62 ± 0.38 






0.70 


0.64- 


0.74 


2.91 ± 0.39 


0.15 


0.13 - 0.17 


0.40 


0.37- 


0.45 


2.30 ± 0.36 






0.50 


0.45 - 


0.54 


2.58 ± 0.32 






0.60 


0.54- 


0.64 


2.35 ± 0.31 






0.70 


0.64 - 


0.74 


2.17 ± 0.36 


0.20 


0.17 - 0.21 


0.26 


0.23 - 


0.30 


2.25 ± 0.33 






0.33 


0.30 - 


0.37 


2.94 ± 0.36 






0.40 


0.37- 


0.45 


2.69 ± 0.33 






0.50 


0.45 - 


0.54 


2.30 ± 0.33 






0.60 


0.54 - 


0.64 


2.81 ± 0.41 






0.70 


0.64- 


0.74 


2.16 ± 0.45 


0.25 


0.21 - 0.27 


0.20 


0.16 - 


0.23 


2.53 ± 0.24 






0.26 


0.23 - 


0.30 


2.46 ± 0.25 






0.33 


0.30 - 


0.37 


2.63 ± 0.28 






0.40 


0.37- 


0.45 


3.04 ± 0.33 






0.50 


0.45 - 


0.54 


2.23 ± 0.31 






0.60 


0.54- 


0.64 


2.13 ± 0.36 


0.30 


0.27 - 0.35 


0.12 


0.08 - 


0.16 


2.55 ± 0.23 






0.20 


0.16 - 


0.23 


2.64 ± 0.24 






0.26 


0.23 - 


0.30 


2.03 ± 0.24 






0.33 


0.30 - 


0.37 


2.24 ± 0.28 






0.40 


0.37- 


0.45 


2.12 ± 0.31 






0.50 


0.45 - 


0.54 


2.61 ± 0.38 


0.40 


0.35 - 0.45 


0.12 


0.08 - 


0.16 


2.73 ± 0.25 






0.20 


0.16 - 


0.23 


2.70 ± 0.28 






0.26 


0.23 - 


0.30 


2.64 ± 0.31 






0.33 


0.30 - 


0.37 


2.53 ± 0.35 






0.40 


0.37- 


0.45 


2.78 ± 0.40 


0.50 


0.45 - 0.58 


0.12 


0.08 - 


0.16 


2.63 ± 0.29 






0.20 


0.16 - 


0.23 


2.87 ± 0.32 






0.26 


0.23 - 


0.30 


2.37 ± 0.34 


0.65 


0.58 - 0.74 


0.20 


0.16 - 


0.23 


1.77 ± 0.34 






0.26 


0.23 - 


0.30 


2.56 ± 0.57 



Acceptance uncertainty: +5/-5 % 



Energy scale uncertainty of ±2 %: -2/+2 % 
Normalization error: 4 % 



Table 17: The fraction of events with a leading neutron at Xl = 0.73 in bins of y and Q 2 for 
the BPC region. 
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0.82 < x L < 1 



Q 2 


Q 2 


y 


y 




ratio (%) 


GeV 2 


range 




ran 


ge 


meas. ± stat 


0.11 


0.10 - 0.13 


0.60 


0.54- 


0.64 


1.15 ± 0.23 






0.70 


0.64- 


0.74 


0.87 ± 0.19 


0.15 


0.13 - 0.17 


0.40 


0.37- 


0.45 


1.19 ± 0.23 






0.50 


0.45 - 


0.54 


1.10 ± 0.19 






0.60 


0.54- 


0.64 


0.82 ± 0.17 






0.70 


0.64 - 


0.74 


0.89 ± 0.21 


0.20 


0.17 - 0.21 


0.26 


0.23 - 


0.30 


1.03 ± 0.20 






0.33 


0.30 - 


0.37 


1.18 ± 0.20 






0.40 


0.37- 


0.45 


1.06 ± 0.19 






0.50 


0.45 - 


0.54 


1.00 ± 0.19 






0.60 


0.54 - 


0.64 


1.16 ± 0.24 






0.70 


0.64- 


0.74 


0.82 ± 0.25 


0.25 


0.21 - 0.27 


0.20 


0.16 - 


0.23 


0.96 ± 0.13 






0.26 


0.23 - 


0.30 


0.79 ± 0.13 






0.33 


0.30 - 


0.37 


1.41 ± 0.19 






0.40 


0.37- 


0.45 


1.03 ± 0.17 






0.50 


0.45 - 


0.54 


0.94 ± 0.18 






0.60 


0.54- 


0.64 


1.09 ± 0.23 


0.30 


0.27 - 0.35 


0.12 


0.08 - 


0.16 


1.47 ± 0.16 






0.20 


0.16 - 


0.23 


1.33 ± 0.15 






0.26 


0.23 - 


0.30 


0.95 ± 0.15 






0.33 


0.30 - 


0.37 


0.93 ± 0.16 






0.40 


0.37- 


0.45 


0.86 ± 0.18 






0.50 


0.45 - 


0.54 


1.14 ± 0.23 


0.40 


0.35 - 0.45 


0.12 


0.08 - 


0.16 


1.33 ± 0.16 






0.20 


0.16 - 


0.23 


1.15 ± 0.16 






0.26 


0.23 - 


0.30 


1.12 ± 0.18 






0.33 


0.30 - 


0.37 


0.93 ± 0.19 






0.40 


0.37- 


0.45 


0.92 ± 0.21 


0.50 


0.45 - 0.58 


0.12 


0.08 - 


0.16 


1.11 ± 0.17 






0.20 


0.16 - 


0.23 


0.94 ± 0.17 






0.26 


0.23 - 


0.30 


1.40 ± 0.24 


0.65 


0.58 - 0.74 


0.20 


0.16 - 


0.23 


0.92 ± 0.22 






0.26 


0.23 - 


0.30 


1.21 ± 0.35 



Acceptance uncertainty: +6/-6 % 



Energy scale uncertainty of ±2 %: +14/-16 % 
Normalization error: 4 % 



Table 18: The fraction of events with a leading neutron at Xl = 0.91 in bins of y and Q 2 for 
the BPC region. 
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e' 




n 

Figure 1: A schematic diagram for the one-pion-exchange model as applied to ep —> e'Xn 
showing the Lorentz invariant variables s', Q 2 , W 2 and t. Also shown is the definition of the 
scattering angle of the positron, 9 e , and the production angle of the neutron, 9 n . In the ZEUS 
coordinate system, the polar angle of the positron is 9 = n — 9 e . 
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Figure 2: A schematic elevation view of the forward proton beam line from the central detector 
to the FNC at Z = +106 m. Magnets are labelled by B (dipoles) or Q (quadrupoles) and a 
number which is their nominal distance from the interaction point in meters. SI through S6 
label the positions of the six leading proton spectrometer (LPS) stations. Between 65 and 80 
m, the proton beam line is bent upwards by 6 mrad by three bending magnets which form a 
momentum analyzer for stations S5 and S6 of the LPS. The neutron exit window is located at 
Z = +96 m, immediately following S6. 
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Figure 3: A schematic view of the FNC, looking towards the interaction region. The calorime- 
ter is segmented vertically into 14 towers. The window of geometric acceptance defined by the 
apertures of the HERA beam-line elements is outlined in towers 7-9. The zero-degree point 
is marked by an x. Towers 11 and 12 have a 10 x 10 cm 2 hole to accommodate the proton 
beampipe. 
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ZEUS 




X (cm) 

Figure 4: Contours of constant neutron production angle compared to the window of the 
geometric acceptance defined by the apertures of the HERA beam-line elements (thick irregular 
curve). The zero-degree point is marked by an x. 
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Figure 5: The kinematic regions in (a) t and (b) p% covered by the angular acceptance of the 
FNC (8 <0.8 mrad) are shown as shaded bands. The solid lines indicate the average t and p 2 T 
values, uncorrected for acceptance, as a function of x L . 
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Figure 6: The slope, b, of the exponential p 
data set. The solid line is given by b(x L ) = 
result obtained using the effective flux of Eq. 
uncertainties only. 



- distribution function of x L for the 1995 DIS 
= (16.32C £ - 4.25) GeV~ 2 . The dashed line is the 
M). The vertical error bars display the statistical 
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Figure 7: (a) Energy spectra, dr/dx L , of neutrons with 9 n < 0.8 mr in PHP, in the intermediate- 
Q 2 region and in DIS. The PHP data have a trigger-correction uncertainty of 2.5% which can 
move them relative to the intermediate-Q 2 and DIS data and, for x L > 0.52, have an additional 
normalisation uncertainty of ±4%. (b) Energy spectra, dr/dx L , of neutrons with 9 n < 0.8 mr in 
DIS for three Q 2 selections. The points in (b) are offset slightly in x L for improved visibility. The 
shaded bands show the systematic uncertainty due to the acceptance of the FNC. There is also 
a normalisation uncertainty of ±5 (±4)% for the PHP (intermediate-Q 2 and DIS) data. The 
energy-scale uncertainty of ±2% is equivalent to a stretch or compression along the abscissa. 
Not shown are the correlated normalisation uncertainties due to the energy-scale uncertainty 
given in Table [I]. 
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Figure 8: (a) The differential cross section for photoproduction of leading neutrons with 9 n < 
0.8 mrad, da LN /dx L , at (W) = 207 GeV. The solid histo gram shows the result of a fit to the 
data using the OPE effective flux factor, / eff of Eq. shown as the dashed histogram, plus a 
background term oc (1 — x L ), shown as the dotted histogram, (b) The normalised cross section 
for 7£> — > Xn at p 2 T = (solid points). Also shown (open circles) are data from the ISR for 
pp — > Xn ||28|| . The curve is the expectation of the one-pion-exchange model for hadron-hadron 
data scaled by 0.41. In both (a) and (b), the statistical uncertainties are shown by the inner 
error bars and the statistical and systematic uncertainties added in quadrature by the outer 
bars. Points for x L < 0.5 (x L > 0.5) are dominated by statistical (systematic) uncertainties (see 
Table §. 
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Figure 9: The ratio r for 0.64 < x L < 0.82 as a function of Q 2 . The ratio r is plotted at the 
average value of Q 2 for the bin. 
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Figure 10: The ratio of the neutron energy spectrum in PHP to that in DIS. The curves 
labelled NSZ [^] and DP [36j are theoretical predictions of the effects of absorption, (Q 2 < 0.02 



GeV )/(Q > 10 GeV ). The vertical bars display the statistical uncertainties. 
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Figure 11: Neutron production (8 n < 0.8 mrad) as a fraction of the inclusive cross section 
in the PHP (Q 2 < 0.02 GeV 2 ), intermediate-Q 2 (0.1 < Q 2 < 0.74 GeV 2 ) and DIS (Q 2 > 4 
GeV 2 ) kinematic regions as a function of Q 2 in bins of y and x L . The PHP point is plotted 
at the average value of (Q 2 ) = 2 • 10~ 3 GeV 2 for the sample. The DIS points are offset in Q 2 
for visibility. The horizontal dotted lines are plotted at the mean r of the plot to guide the 
eye. Only statistical uncertainties are plotted for the DIS and intermediate-Q 2 regions. For the 
PHP point, the uncertainty plotted is the statistical uncertainty added in quadrature with the 
PHP relative normalisation uncertainty. 
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Figure 12: Neutron production (6 n < 0.8 mrad) for the intermediate-Q 2 region, 0.1 < Q 2 < 
0.74 GeV 2 , as a fraction of the inclusive cross section and as a function of y for the low 
(0.2 < x L < 0.64), medium (0.64 < x L < 0.82), and high (0.82 < x L < 1) x L ranges. The 
dotted lines show the mean values of the ratio for each x L range. Not shown are the correlated 
systematic uncertainties given in Table |I[ 
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Figure 13: Neutron production (9 n < 0.8 mrad) for the DIS region, Q 2 > 4 GeV 2 , as a fraction 
of the inclusive cross section and as a function of x for the low (0.2 < x L < 0.64), medium 
(0.64 < x L < 0.82), and high (0.82 < x L < 1) x L ranges, in the indicated bins of Q 2 . The dotted 
lines show the result of fitting a power law in x to the ratio. The solid lines show the result 
of a fit to the ratio linear in both lnx and InQ 2 , as discussed in the text. Not shown are the 
correlated systematic uncertainties given in Table 0. 
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Figure 14: F 2 LN(3) for 6 n < 0.8 mrad for 0.20 <x L < 0.64, 0.64 <x L < 0.82 and 0.82 <x L < 1.0 
in the intermediate-Q 2 region, 0.1 < Q 2 < 0.74 GeV 2 , as a function of x. The uncertainties 
shown are statistical added in quadrature with the uncertainties from the ZEUS Regge fit to 
F 2 . The curves show F 2 (x, Q 2 ) for the proton scaled by the average value of r divided by the 
bin width Ax L . Not shown is an additional correlated systematic uncertainty of ±8.8% from 
the acceptance, the energy scale and the overall normalisation uncertainty. 
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Figure 15: F 2 w for 6 n < 0.8 mrad for 0.40 <x L < 0.46, 0.70 < x L < 0.76 and 0.88 < x L < 1.0 
in the DIS region, Q 2 > 4 GeV 2 , as a function of x. The uncertainties shown are statistical 
added in quadrature with the uncertainties from the ZEUS NLO QCD fit. The curves show 
F 2 (x,Q 2 ) for the proton scaled by the average value of r divided by the bin width Ax L . Not 
shown is an additional correlated systematic uncertainty of ±8.8% from the acceptance, the 
energy-scale and the overall normalisation uncertainty. 
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Figure 16: A comparison of the ZEUS and HI values of at (a) x L = 0.3, where the 

p T ranges covered by the two measurements coincide, and (b) x L = 0.7, where the ZEUS data 
have been adjusted to the transverse-momentum range covered by HI, p T < 0.2 GeV, using the 
form of Eq. (0). For ZEUS, the statistical uncertainties are added in quadrature with that 
from the ZEUS NLO QCD fit. A normalisation uncertainty of 6.7% (7%) remains at x L = 0.3 
(0.7) due to the acceptance, energy scale and overall normalisation uncertainties. For HI, the 
statistical and systematic uncertainties were added in quadrature. The uncertainties on the HI 
points are dominated by the x L -dependent uncertainties and so are strongly correlated. There 
is an additional normalisation uncertainty of 5.7% (not shown). 
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Figure 17: (a) The ry max distribution in DIS for neutron-tagged events (points). The histogram 
shows the inclusive 7/ max distribution, normalised to the neutron-tagged data for 7/ max > 1.8; (b) 
the fraction, Rlrg> of leading neutron events with Q 2 > 4 GeV 2 passing the LRG criterion 
T] m&K < 1.8, as a function of x L . 
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Figure 18: F£ as a function of for the pion in bins of Q 2 determined for 0.64 < x L < 0.82, 



using Eq. ([0J). The pion flux used to determine F% is the effective OPE flux (EF) used in 
hadron-hadron charge-exchange reactions, Eq. (f|). The uncertainty shown on F^' EF arises from 
the statistical uncertainty due to the leading neutron added in quadrature with the uncertainty 
on i<2. Not shown are the correlated systematic uncertainties given in Table [I]. The dotted 
lines are F 2 (x,Q 2 ) for the proton, scaled by 0.361. The solid curves are F£ from the GRV 
parameterisation 
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Figure 19: F% as a function of for the pion in bins of Q 2 determined for 0.64 < x L < 0.82. 
The pion flux used to determine F% is the flux obtained using the additive quark model (AQM) 
of Eq. (pip . The uncertainty shown on F%' arises from the statistical uncertainty due 
to the leading neutron added in quadrature with the uncertainty on F 2 . Not shown are the 
correlated systematic uncertainties given in Table |I[ The solid curves are F£ from the GRV 
parameterisation [7f| while the dotted curves are from the Sutton et al. parameterisation |77 . 
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